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04 ' A basic issue for Navier-Stokes (NS) fluids is their characterization in terms of the so-called NS 

}_( ' phase-space classical dynamical system, which provides a mathematical model for the description 

C^ , of the dynamics of infinitesimal (or ideal) tracer particles in these fluids. The goal of this paper is 

to analyze the properties of a particular subset of solutions of the NS dynamical system, denoted 
as thermal tracer particles (TTPs), whose states are determined uniquely by the NS fluid fields. 

[~^ \ Applications concerning both deterministic and stochastic NS fluids are pointed out. In particular, 

in both cases it is shown that in terms of the ensemble of TTPs a statistical description of NS fluids 
can be formulated. In the case of stochastic fluids this feature permits to uniquely establish the 
corresponding Langevin and Fokker-Planck dynamics. Finally, the relationship with the customary 
statistical treatment of hydrodynamic turbulence (HT) is analyzed and a solution to the closure 
problem for the statistical description of HT is proposed. 
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CO . A fundamental aspect of theoretical fluid dynamics is represented by the discovery of the thermal tracer particles 

•^ . (TTPs) for Navier-Stokes (NS) fluids recently reported (see Ref.[l|)- The latter represent a suitable subset of the 
f— ^ ,, so-called ideal tracer particles (ITPs [2]) and are defined in such a way that their states are uniquely dependent, in 
PsJ ■ a sense to be specified below, on the local state of the fluid. A basic implication of the result is that an appropriate 
T-H ' statistical ensemble of TTPs should reproduce exactly the dynamics of the fluid. In other words, it should be possible 
^ to determine the fluid fields characterizing the fluid state by means of suitable statistical averages on the ensemble of 
TTPs, and in particular performed so that they satisfy identically a required set of fluid equations. The conclusion 
is expected to apply, in principle, to arbitrary NS fluids described as mesoscopic, i.e., continuous fluids, which can be 
either viscous or inviscid, compressible or incompressible, thermal or isothermal, isentropic or non-isentropic. 

We shall assume, for this purpose, that the state of these fluids is represented by an ensemble of observables {Z{y, t)} 
= {Zi{Y,t),i = 1, ..,n} (with n an integer > 1), i.e., fluid fields, which can be unambiguously prescribed as continuous 
and suitably smooth functions, respectively, in fi x / and in the open set fi x /, with fi C M^ and / = M being the 
configuration space and time axis respectively. We intend to show that, as a basic consequence, the Newtonian state of 
each TTP, namely x(i) = x = (r, v), with r and v denoting respectively the particle position and velocity, is advanced 
in time in terms of a suitable acceleration field F = F(x, i) which can be defined in such a way to depend only on 
the state of the same particle {mean- field acceleration). Remarkably, it is found that F can be uniquely prescribed 
in such a way to determine self-consistently the time evolution of the complete set of fluid equations characterizing 
the fluid. This implies that TTPs must reproduce exactly the dynamics of the fluid. In other words, by means of 
appropriate statistical averages on the ensemble of the TTPs, it is possible to uniquely determine the time-evolution 
of the fluid state, in such a way that it satisfies identically the required set of fluid equations. 



A. Lagrangian dynamics of ideal tracer particles 

A key aspect of fluid dynamics is tfie proper definition of the phase-space Lagrangian dynamics for continuous fluid 
systems, whereby possibly all the fluid fields characterizing the actual fluid state {Z{r,t)} can be identified with 
suitable statistical averages on appropriate ensembles of (fictitious) particles. Thus, for example, in the case of an 
incompressible NS fluid, this would require to represent both the fluid velocity Y{r,t) and the fluid pressure p{r,t) 
in terms of suitable statistical averages of an appropriate probability density. This goal can be realized by means of 
the inverse kinetic theory (IKT) developed by Ellero and Tessarotto (see Refs.^3H5|]). This refers, in particular, to 
the phase-space dynamics of ideal tracer particles, namely rigid extended classical particles immersed in the fluid, 
all having the same support and infinitesimal size such that during their motion they do not mutually interact and 
do not perturb the state of the fluid. Depending on their incrtial mass mp, ITPs can belong to different species 
of particles; thus, in general, their mass can differ from that of the corresponding displaced fluid element mp- On 
the other hand, ITPs carrying the mass mp = mp will be denoted as the NS ideal tracer particles (NS-ITPs). In 
the following, in order to characterize the Lagrangian dynamics of NS fluids, ITPs will be identified only with NS- 
ITPs. In this framework, it follows that ITPs can undergo, by assumption, only "unary" interactions with external 
force-fields and with the continuum fluid. Namely, in both cases they are subject only to the action of a continuum 
mean-field acceleration which depends only on the local state of each particle. As a consequence, ITPs can be treated 
as Newtonian point-like particles characterized by a Newtonian state x = (r, v) spanning the phase-space T = il x U, 
with the position r and the kinetic velocity v belonging respectively to the configuration space of the fluid J7 (in the 
following to be identifled with a bounded subset of R'^) and the velocity space [/ = R'^. 

B. The Navier-Stokes dynamical system 

By assumption, the state x of a generic ITP advances in time by means of a Newtonian classical dynamical system 

(DS) deflned in terms of the vector fleld X(x, t) = <v, ;;^K = F > , with ;;^K = F a suitable mean-fleld acceleration. 

This is identifled with the flow (Tt^t) generated by the initial value problem associated to the deterministic equations 
of motion (Newton's equations) 

i^ = X(x,i), 

x(io)=Xo. ^'> 

Such flow is referred to as Navier-Stokes dynamical system (NS-DS) and is a homeomorphism in P with existence 
domain P x /, of the type 

Tt^,t ■ Xo -^ x(i) = Tt^^yio, (2) 

with t G / C M and Tt^.t being a measure-preserving evolution operator associated to X(x, t). Thus, by definition, the 
NS-DS is uniquely prescribed by the couple {x, X(x, t)}, with x(t) = x = (r, v) to be identified with the instantaneous 
Newtonian state of a generic ITP. 

C. The relative-dynamics NS-DS 

The state of a TTP can be equivalently represented in terms of the relative-dynamics Newtonian state y — (r, u), 
with u(i) = v(i)— V(r(i), t) denoting the relative kinetic velocity defined with respect to the local fluid velocity V(r, t). 
As a consequence, by introducing the local phase-space diffeomorphism x = (r, v) — >■ y = (r, u), the dynamical system 
([2]) can be cast in terms of y = (r, u). Then, © can be equivalently represented by the homeomorphism in P: 

TiT--yo^yit)=T,'l'f'yo, (3) 

to be identified with the relative-dynamics NS dynamical system (RD-NS-DS), T^ ^ being the flow generated by 
the initial-value problem 

" ^ = u + V(r,i), 

f=F4x,i), (4) 

y[to) = Yo. 



Here, for a prescribed form of the mean-field F(x, i) (see below), Ftj(x, i) denotes the kinetic relative acceleration 
which is defined as 

F„(x, t) = F(x, t) - FH(r, i) - u • VV(r, i), (5) 

with F/f (r, t) being the NS fluid acceleration defined by Eq. (|142p (see Appendix A). To establish on rigorous grounds 
the connection with the corresponding fluid description, the vector field F(x,i) must be suitably determined. For 
definiteness, we shall consider here the case of a Navier-Stokes thermofluid described either by the compressible or 
incompressible Navier-Stokes- Fourier equations, requiring that the fluid fields {Z(r, t)} are strong solutions either 
of the compressible or incompressible Navier-Stokes-Fourier problems (CNSFE or INSFE problems respectively; see 
Appendix A). In the case of CNSFE the state of the fluid is defined by the ensemble of smooth fluid fields 

{Z(r, t)} = {p(r, i), V(r, i),p(r, t), T(r, t), Srit)} , (6) 

with p(v,t) > 0, V(r, i), p{r,t) > 0,r(r,i) > and Sxit) denoting the fluid mass density, the fluid velocity, the fluid 
scalar pressure, the fluid temperature and the global thermodynamic entropy respectively. In particular, by introducing 
an arbitrary reference mass m > 0, for example to be identified with the particle mass mp defined above, the notion 
oi fluid number density n{r,t) = —p{r,t) can be introduced. As an alternative, in the following the set {Z(r,t)} will 
be replaced by the reduced set of fluid fields 

{Z,{r,t)} = {p{r,t),V{r,t),p,{r,t),ST{t)} . (7) 

Here pi(r,t) denotes the kinetic pressure, i.e., a strictly positive scalar observable defined as 

Pi(r, t) - po(t) + p{r, t) ~ 0(r, t) + n{r, t)T{r, t), (8) 

with po(i) being the pseudo-pressure and 0(r, t) the potential associated to the conservative part of the volume force 
(see Appendix A) [55| . For later use we introduce here also the notion of specific kinetic pressure pi{r,t): 

pi{r,t)=pi{r,t)/p{r,t). (9) 

In particular, we shall assume that in Jl x / both (j>{r,t), po{t) and Srit) are uniquely defined. In particular they are 
such that: 

• <j>{r,t) is bounded and such that, for r = Rq 

(/.(Ro,i)|-0, (10) 

with Ro being a suitable vector belonging to fl. 

• po{t) is a suitably-prescribed smooth real function defined so that pi(r, t) remains strictly positive in O x / (see 
Axiom ^2, Section 4). 

• The initial value po(to) is in principle an arbitrary constant to be prescribed in such a way that po (tg) > poinf, 
with poinf being such that pi{r,to) vanishes locally in H, (see Axiom ^2, Section 4). 

D. The IKT statistical description 

It must be stressed that the vectors v and u defined above can be interpreted as stochastic variables [see Appendix 
B] , and the corresponding equations of motion, Eqs. ([T]) and Q viewed as Langevin (i.e. , stochastic) equations, provided 
a statistical description in terms of a suitable probability density is introduced for them. Statistical descriptions of 
this type, based on classical statistical mechanics (CSM), can be adopted in principle both for classical and quantum 
fluids (see for example Refs.fj-Q and Q) characterized either by deterministic or stochastic flows [3, [i3l- This is 
realized by introducing an appropriate axiomatic approach denoted as statistical model, represented by a set {/, F} , 
with f{'x.,t) denoting a suitable kinetic distribution function (KDF) - or a probability density function (PDF), to be 
identified with the so-called 1-point PDF - which is defined in the phase-space F. Its is worth noting that such a type 
of approach can be determined in accordance with the GENERIC dynamical model {/, F} developed by Grmela and 
Ottinger [ll|, II4I. This is defined in such a way to prescribe: 

• A bundle structure on F [ll|, i.e., a mapping between F itself and the ensemble of fluid fields {Z} (or {Zi}), 
generated via appropriate statistical averages, i.e., phase-space moments, of /(x, t). 



• A functional class {/} to which /(x,i) belongs M^, ^. This should be based exclusively on the knowledge of the 
same set of fluid fields. The prescription of {/} includes also the initial condition on /(x, i) at i = to- 

• A suitable phase-space dynamics: this is introduced in terms of the NS-DS [see Eq.(l2|)]. This should inherit 
the basic properties of the fluid system, i.e., in particular, the conservation of mass and momentum, the energy 
balance equation and the entropy law. Thus, in the case it is described by deterministic, dissipative and 
irreversible fluid equations, it should be a deterministic, non-conservative and irreversible dynamical system. 

The problem of its construction (Frisch, 1995 [i3)j i-e., the actual definition of the vector field F, has remained for 
long time unsolved (see for example, Vishik and Fursikov, 1988 14 [ and Ruelle, 1989 [i^|, where approaches based 
on kinetic theory were attempted). Nevertheless, in the pas t various models for the dynamics of tracer particles 
in incompressible fluids have actually been developed |16l420| . which should manifestly apply, at least in principle, 
also to ITPs. As a consequence, their dynamics should be controlled only by the mean-field force F produced by 
the unperturbed fluid flelds. Since the original Basset-Boussinesq-Oseen approach |2l|-(23j formulated in the case 
of a uniform flow, several attempts to evaluate the form of F, as well as the vector field F„, have appeared |16l - 
[201 • AH such approaches propose "ad hoc" modifications or corrections of the same equation, exclusively based on 
phenomenological arguments, in order to adapt it for the treatment of non-uniform flows in NS fluids. A popular 
form (for F) frequently adopted for fluid simulations is the one developed by Maxey and Riley (1982 [20])- These 
treatments appear questionable because of the critical common assumption on which they are based. Precisely, the 
requirement that the tracer-particle velocity v(i) (kinetic velocity) remains always suitably close to the local velocity 
of the fluid V(r,t) evaluated at the position of the moving particle r = r(t). This implies that for all ITPs the 
asymptotic condition 



|u(i)|«|V(r(t),i)l 



fll) 



should hold at any time t. As a consequence, in the determination of F all contributions proportional to the relative 
velocity should be considered negligible. The constraint (jlip imposes, however, potentially serious limitations on 
particle dynamics. In fact, it may easily be violated either due to the arbitrariness of the particle initial conditions 
[in Eq.([T])] - in fact ITPs can be injected in a fluid with arbitrary initial velocity - or because ITPs initially at rest 
with respect to the fluid might develop finite relative velocities, thus causing Eq. pTj) to fail. On the other hand, 
the determination of the exact ITP dynamics is manifestly of fundamental importance in order to obtain detailed 
quantitative theoretical and numerical predictions in fluid dynamics. 

A first-principle solution of this problem has recently been proposed in Ref. [2| adopting a suitable statistical de- 
scription for incompressible Navier-Stokes thermofluids (see Appendix A), i.e., the representation of the dynamical 
system ^ in terms of the IKT- statistical model {/(x, t), F}, obtained in the framework of the so-called inverse kinetic 
theory (IKT; see Refs.tJ,Q). In such a formulation. 



/(t) = /(x,t) 



(12) 



is identifled with a KDF whose velocity and phase-space moments are prescribed in terms of a suitable subset of the 
fluid fields {Z}, to be identified with the ensemble {Zi}. Hence, provided J dvf{x,t) > 0, the corresponding velocity 

u 
PDF is 

J[x,t) = f{x,t)/ j d^f[x,t). (13) 

u 
The KDF / is required to satisfy in F x / the statistical equation 

Lf = 0, (14) 

denoted as inverse kinetic equation (IKE [5,]) in Eulerian form, with L = -^f + v • V/ + ^ • (F/) being the Liouville 
streaming operator. The same equation can be equivalently cast in terms of the integral Lagrangian IKE 



/(x,t) = L(r,,t„x) 



dTt,t^^ 



dy 



Here foi^o) is a suitable initial KDF, while 

5T.t„x 



9x 



(. to 



(15) 



(16) 



(Liouville theorem). In addition, the following assumptions are introduced (see points a,b,c,d below): 

a) A particular solution of Eg. (1141) is provided by the local Gaussian distribution function (kinetic equilibrium) 

/M(x,t) = ^^exp{-X2}, (17) 

(tt) 2 Vff^ 

where X'^ = :^jfj^ and Vth = y ~T" denotes the thermal velocity. 

b) In terms of a suitable set of the velocity and phase-space moments, the complete set of fluid equations are 
constructed from IKE. 

c) Let us assume that /(x, t) is strictly positive in F and admits for all (x, i) G F x / the Boltzmann- Shannon (BS) 
statistical entropy |S] (also known as differential entropy). This is defined as the functional 

SU{t)) = -al I dx/(x,i) In /(x,i) + c,, (18) 

r 

where ai ^ and Ci are arbitrary real constants independent of (x, t), to be suitably defined. In information theory 
the Boltzmann-Shannon entropy S{f(t)) can be intended as a measure of the ignorance on f{t). Here we remark 
that, by suitable definition of the constants ai and ci, S{f{t)) can be represented in terms of the corresponding BS 
entropy for the phase-space PDF f{t) = f{t)/AI, where M — J dxf{ii.,t). In fact, letting 

r 

S{f{t)) = -J dx7(x,t) ln7(x, t), (19) 



5(/(i)) = ^5(/(i))+lnM, (20) 



it follows that 



which is again of the form (TTSl) . In the following we shall set in particular ai = 1 and ci = in Eq. lfTSl) . 

d) The time-derivative of the pseudopressure dpo{t)/dt is determined by suitably prescribing the entropy production 
rate f 5(/(i)) M- 

As a basic consequence of the previous assumptions it is follows that [5|, Q : 

1. The mean-field F is generally functionally dependent on /, i.e., it is of the form F = F(x, t; /). 

2. F is defined up to an arbitrary real gauge vector-field AF obeying the gauge condition 

^•[AF(x,i;/)/(x,i)]=0. (21) 

3. The choice of the gauge field AF does not affect the time- evolution of f{x,t). Consequently, the corresponding 
velocity and phase-space moments of the IKE are in all cases necessarily unique. 

4. In the case / = /m and up to the gauge field AF, the functional form of the mean field F(x, t; fm) is uniquely 
determined. 

5. In the case / 7^ f^ the determination of F, again up to the gauge field AF, requires suitable kinetic closure 
conditions. In fact, in principle, in such a case F might depend on arbitrary higher-order phase-space moments 
of / which vanish in the case / = /a/. To this end, in the case of incompressible NS fluids, in Refs.0, M, @| it 
was assumed that F(x, t; f) (and AF) can be represented as polynomials of lowest possible degree with respect 
to the relative kinetic velocity u and depend on the lowest-order and minimal number of velocity moments of 
the KDF. 

We remark that, although in the context of IKT the unique specification of AF (and hence of F) is superfluous, 
its determination is, instead, manifestly required in order to uniquely prescribe the dynamics of ITPs. In particular 
in Ref. f3| , based on the analogy with extended thermodynamics [6] , AF was identified with a first-degree polynomial 
with respect to u of the form AF — | VV • u — ^u • VV. In the following we intend to propose a generalization of IKT 
for compressible Navier-Stokes thermofluids (Section 2) and in which AF is uniquely determined, based on suitable 
physical assumptions. For this purpose, leaving initially unspecified the form of the gauge field AF, we intend to 
prove that AF is uniquely prescribed imposing the requirements stemming from the following Gedanken (conceptual) 
experiment (GDE). 



II. GEDANKEN EXPERIMENT 

For a prescribed continuous fluid system, such as a compressible/incompressible NS thcrmofluid, the problem arises 
whether there might exist a subset of the ensemble of ideal tracer particles (ITPs) for the dynamical system ([2]) 
such that their Newtonian state and corresponding time evolution depend only on the state of the fluid {Z}. In the 
following the subset of ITPs which exhibit these properties are referred to as thermal tracer particles (TTPs) . Such 
a result was reached for incompressible and isothermal NS fluids in Rcf. T| . Here we claim that it should be possible 
to extend the same conclusion to arbitrary compressible and non-isothermal NS fluids by performing a conceptual 
experiment {Gedanken experiment) on such a type of fluid, i.e., looking at the properties of the IKT-statistical models 
{/, F}. The conjecture is suggested by the following arguments: 

• The state of the fluid is solely dependent on the fluid flelds, which in the case of a compressible NS thermofluid 
can be identified with the set {Zi}. 

• The time-evolution of {Zi} as determined by CNSFE is necessarily independent of the KDF f{x.,t) and of the 
NS-DS ([2|). In fact, obviously the CNSFE (or INSE) problem cannot depend on the functional form of f{x,t). 

• On the other hand, in the context of IKT, the time evolution of the KDF is determined by the Liouville 
operator L, which enters the corresponding Liouville equation. As pointed out in Refs.0, la], this generally 
contains a vector field F whose form can depend functionally also on the same KDF /(x, t), namely is of the 
form F = F (x, i; /). Despite this, the time evolution of the fiuid fields {Zi} generated in terms of the KDF by 
the Lagrangian IKE (fT5|) through the NS-DS ([2]) , can easily be shown to remain independent of the functional 
form of the same KDF. 

A. GDE requirements 

On the basis of these considerations, here we conjecture that TTPs should exist as a subset of ITPs and fulfill the 
following properties: 

1. GDE-requirement #1: their time evolution, as determined by the vector field F in terms of the NS-DS ©, 
should remain at all times t G / independent of the particular form of the KDF /(x, i). As a consequence, 
for them the form of the mean-field force F should be also independent of the KDF f{x,t) [introduced in the 
IKT-statistical model {/, F}], namely simply of the form F = F (x, t). 

2. GDE-requirement #2: for prescribed initial conditions, their Newtonian states x(t) = x = (r, v), and equiva- 
lently also y{t) = y = (r, u), should depend solely on the fluid fields {Zi}. 

In addition, one should expect that for all TTPs: 

3. GDE-requirement #3 - Local magnitude of u(i): the magnitude of their instantaneous relative velocity |u(i)| = 
lu(r,t)| remains at all times t e / proportional to the local thermal velocity wt/i(r,i), i.e., of the form 

\u{t)\^^vth{r,t), (22) 

with pi(r, t) > 0, r = r{t) and (3 denoting respectively the kinetic pressure ([5]), the instantaneous position of the 
same particle and an appropriate non- vanishing constant, i.e., a function independent of (r, t). This means that 
/3 is necessarily determined by the TTP initial condition (see discussion below, after TIIM.2); 

4. GDE-requirement #4 " Kinetic constraint on the local direction of u(i): let us introduce for u{t) the represen- 
tation 

u{t) ^ Pvth{v,t)n{r,t), (23) 

with n(r, t) being the unit vector prescribing the local direction of u(t). Then, if at a given point the constraint 
([22)l is satisfied, in order to warrant that u(<) satisfies it also at time t -\- dt (with dt being infinitesimal), it 
is necessary to require that the unit vector n(r,t) be tangent to the local isobaric surface pi{r,t) — const. 
Therefore, for a non-uniform kinetic pressure satisfying locally Vpi(r, t) ^ 0, the unit vector n(r, t) must satisfy 
the kinetic constraint: 

n(r,i)-Vpi(r,t)=0. (24) 



As a consequence, the direction of u(<) is necessarily uniquely determined, once the initial conditions ([T]) and 
consequently its initial direction 

n(r,io) = no(r) (25) 

have been set. Hence, for a non-uniform specific kinetic pressure pi, the unit vector n(r, i) must be orthogonal 
to the unit vector 

, / ,N Vpi(r,t) 

i.e., the kinetic constraint 

n(r,t)-b(r,t)=0 (27) 

must hold identically for all (r, t) £ fi x /. 

5. GDE-requirement #5 - Time evolution of n(r,t): the unit vector n{r,t) satisfies an initial- value problem of the 
form 

^^=nir,t)xnir,t), ^28) 

n(r(t„),io) = n(ro,to), ^ ^ 

with ft{r,t) denoting a suitable pseudo- vector. Without loss of generality we shall require that fl{r,t) is a 
smooth real vector function defined in fi x / and that it is defined also in the limit pi{r, t) ^> 0+. 

6. GDE-requirement #6 - Rotation dynamics of n(r,i): we require that the unit vector n(r, i) exhibits a rotation 
motion with respect to the direction b(r,t) which is determined by the parallel component of the vorticity. In 
other words, the n{r,t) is required to satisfy the constraint 

n(r,t)-b(r,i) = -e(r,t).b(r,i), (29) 

where ^(r,t) = V x V(r, i) is the vorticity field. As a consequence of Eqs. ([28l) and ([29|). the particle relative 
velocity exhibits a rotation caused by two distinct physical mechanisms. The first one is due to the rotation of 
the unit vector b(r, t) characterizing the isobaric surfaces, while the second one is due to the intrinsic rotation 
of the parallel component of fl{r, t) around b(r, t) as determined by the local vorticity field. In particular, the 
constraint placed by Eq. (j29p establishes a direct connection between TTP dynamics and fiuid vorticity and has 
important implications for the treatment of strong turbulence in the framework of TTP statistics (see related 
discussion in subsection VII. B below). 

III. GOALS OF THE INVESTIGATION 

In this paper we point out that IKT can be determined in such a way to satisfy the requirements dictated by the 
GDE. For this purpose, first the IKT statistical description {/, F} earlier pointed out in Refs. 0, S la] is extended 
to the treatment of compressible and non-isothermal NS fluids. Next, the NS dynamical system is shown to admit 
particular solutions which arc of the form of TTPs, namely ITPs for which the particle state takes the form indicated 
above [see Egs. d^ and (P71) - (P5)) ]. More precisely, extending the results pointed out in Ref.[l[, and holding in the 
case of incompressible and isothermal NS fluids, here we intend to prove that for compressible, non-isothermal fluids 
satisfying the CNSFE Problem (see THM.l in Section 4): 

• Goal jj=l - TTPs are particular solutions of the NS dynamical system (see Section 5, TIIM.2). 
In such a setting, we claim that the following additional properties are fulfilled: 

• Goal #2 - For all TTPs a unique realization exists for F -|- AF satisfying the requirements of GDE. 

• Goal #3 - In terms of the ensemble of TTPs a reduced- dimensional statistical model {/i,Fi} {TTP -statistical 
model) is introduced (Section 6, TIIM.3). 



Another interesting application concerns the treatment of stochastic flows. For this purpose the fluid fields p6)) are 
assumed to admit, in terms of suitable stochastic variables a — {ai,i = l,k} G Va C R*^, with fc > 1, a stochastic 
representation of the form [2, [ifll 

{Zi} = {Zi(r,t,a)}, (30) 

to be defined in terms of a suitable stochastic model {g{r,t,a),Va} (see Appendix B). We intend to show that for 
incompressible fiuids: 

• Goal #4 - The Langevin equations associated to TTPs dynamics provides a possible mathematical model for 
tracer-particle motion in the presence of fluctuating fluid fields (Section 7, subsection 7.1). 

• Goal ^5 - The stochastic-averaged KDF of the TTP- statistical model {/i, Fi} satisfies a Fokker-Planck statistical 
equation and a H-theorem (Section 7, subsections 7.2 and 7.3). 

• Goal =1^6 - The IKT and TTP statistical models {/, F} and {/i,Fi} are suitably related to the customary 
statistical treatment of turbulence due to Hopf Rosen and Edwards (HRE approach [24-r26l |') (see Section 7, 
subsection 7.4). 

• Goal 7^7 - The TTP-statistical model {/i, Fi} provides a solution to Closure Problem for the statistical descrip- 
tion of hydrodynamic turbulence (HT) (Section 7, subsection 7.4). 

IV. IKT FOR COMPRESSIBLE NS THERMOFLUIDS 
A. Axiomatic formulation 

The basic requirements of the IKT statistical model {/(x, i),F} have been discussed elsewhere @, Q- In the case of 
a compressible NS thermofiuids, these can be re-formulated as follows. First we require that the KDF f{t), solution of 
IKE [see Eq. (fT4| ]. uniquely determines the complete set of fiuid fields {Zi}, in terms of suitable phase-space moments 
of the same KDF. This is obtained imposing the following axiom: 

Axiom #1 - Correspondence principle: For compressible (or incompressible) NS thermofiuids - in the closure of the 
fluid domain fi, where by definition p(r, i) > in fi x / - the following functional constraints hold 

/dv/(x,i)=p(r,t), 
u 
-r^ fdwwf{^,t) = V(r,t), 

/dvi7/2/(x,i)=Pi(r,t), 

S{f{t)) = Srit), 

which are referred to as correspondence principle for {Zi}. Here we remark that due to the arbitrariness of the 
constants ai and ci appearing in the definition of the BS entropy [see Eq. ()18|) ] the last equation can also be replaced 
by S{f{t)) = a|5T(i) + C2, with a2 ^ and C2 being two arbitrary real constants independent of (x,i). Furthermore, 
we shall require that f{t) admits also the higher-order moments 

Q(r,i) = Jdvu^f{r,u,t), .„„^ 

n(r,i) = /dvuu/(r,u,i), ^^^^ 

to be denoted as extended fluid fields. As a consequence, by an appropriate definition of the mean-field, the corre- 
spondence principle and IKE [see Eq. (fT4|) ] must deliver the complete set of fiuid equations, i.e., respectively CNSFE 
or INSFE {fluid closure condition). 

Second, consistent with GSM and the second principle of thermodynamics, the initial KDF f{to) = /o(x) and the 
pseudopressure po (i) are uniquely prescribed. This is obtained introducing the second axiom: 

Axiom #2 - Entropic principle: 

This consists in the following three requirements: 

A) At the initial time to € I the initial KDF f{to) = fo{^) is determined in such a way to maximize the BS-entropy 
S{f{to)) in a suitable functional class {/(to)}- This coincides with the axiom of GSM known as principle of entropy 
maximization (PEM, Jaynes 1957 |27j). 



B) The time derivative of pseudo-pressure po (t) is prescribed for all t G I in such a way that the entropy law (|13ip 
[see Appendix A] is identically fulfilled. 

C) The initial condition po{to) is determined by suitably prescribing the initial value of the BS entropy S{f{to))- 
The first requirement is met as follows. Denoting by S the Frcchet functional derivative operator, let us introduce 

the first and second variations oi S{f{to)), SS{f{to)) and 6'^S{f{to)). Then, the initial KDF fo G {/(^o)} is determined 
imposing the variational equation 

SSifito))\^^^0 (33) 

subject to the inequality 

S'S{f{to))\^,^<0. (34) 

The determination of dpo{t)/dt is obtained, instead, in such a way to warrant - for consistency with the correspondence 
principle and the entropy law - that the weak H-theorem 

^5(/(i)) = I^T(i) > (35) 

holds for all t £ I. In the case of isentropic fiows this reduces to the constant H-theorem: 

|5(/(i))-|5T(t)=0. (36) 

Furthermore, the initial kinetic pressure pa (to) is prescribed in such a way that the initial BS entropy S{f{to)) 
vanishes, i.e., 

Sifito)) = 0. (37) 

The initial condition (j37p on the BS entropy is equivalent to demand that the measure of ignorance S{f{to)) is zero. 
Since S{f{to)) < S{fM{to)) this means that the Gaussian KDF fuito) must admit the BS entropy and hence that 
the kinetic pressure is necessarily pi{r,to) > 0. Due to the arbitrariness of po (to) this requirement can always be 
satisfied. 

Finally, suitable kinetic closure conditions are introduced to determine F(x, t; f) in the non-Gaussian case [/(x, t) j^ 

/M(x,i)]: 

Axiom #3 - Kinetic closure conditions: For this purpose, in analogy with INSE |5|, we shall assume that F(r, u, t; f) 
is a polynomial of lowest possible degree with respect to the relative kinetic velocity u and depends on the lowest-order 
and minimal number of velocity moments of the KDF. In particular: 

- Axiom #3a) F(r,u,t; f) can depend, besides {Zi}, only on the minimal set of extended fluid fields {Q,n}; 

- Axiom #3b) F(r, u, t; f) depends only linearly with respect to Q and H. 

Let us briefly comment on these requirements. In principle the vector field F(r, u, t; /) might depend on arbitrary 
higher-order moments of the KDF. In fact, due to the arbitrariness in its definition (see discussion above), it is 
always possible to include an additive contribution which vanishes identically in the case / = /m (Gaussian KDF) 
and which does not contribute the velocity moments of the Liouville equations corresponding to the weight functions 
G = (l,v, u^/3). Hence the above closure conditions warrant that such contributions are excluded, so that Axioms 
#3a and #3b actually realize the minimal requirements when / is non-Gaussian. 

The motivation for the precise choice of Axiom #3 is mathematical simplicity. In fact, within the framework of 
IKT, the mean-field F is not a physical observable, both for Gaussian and non-Gaussian KDFs, and therefore it 
remains intrinsically non-unique. Its indeterminacy in the case of a non-Gaussian KDF arises because the vector 
field F(r, u, i;/) might include in principle higher-order velocity moments of the KDF, besides the fluid flelds {Zi} 
which are by construction the only observables. These additional moments have no physical meaning (i.e., they are 
not observables), and therefore remain completely undetermined in the framework of IKT. Hence, r(r, u, t; f) should 
depend only on a minimum finite number of fiuid fields which are required for the validity of the theory. This means 
that there must exist a finite subset of moment equations which coincide with CNSFE. To further clarify the issue, 
we notice that F(r, u, t; /) can always be given a polynomial representation in terms of the relative velocity u. Such 
terms would necessarily depend on higher-order velocity moments which vanish in the case of the Gaussian KDF and 
can always be prescribed in such a way not to contribute to the same moment equations. Unless Axioms #3a and 
#3b are introduced, such additional contributions to r(r, u,t;/), depending, besides the CNSFE fluid flelds, also on 
Q and n, would remain undetermined. Due to the intrinsic freedom of their choice, in the following they will be 
set identically equal to zero. As clarifled below, such an assumption is equivalent to require Axioms #3a and #3b. 
This choice does not affect the validity of the CNSFE problem and does not constraint in any way its solutions. In 
conclusion, in view of these considerations. Axiom ^3 can be viewed as a set of kinetic closure condition which are 
needed for the prescription of the mean- field F(r, u, t; /) and the related kinetic equation (ITi|) . 
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B. The IKT statistical model for compressible thermofluids 



Based on the axiomatic formulation given above [Axioms ^1-^3], we can now proceed to the expHcit determination 
of the mean-field F = F(x, i, /) appropriate for a compressible thermofluid satisfying the CNSFE problem [defined 
by Eqs. p28p - (|13ip and the initial-boundary conditions (I163P : see Appendix A]. It is immediate to show that in such 
a case the form of the mean-field F(x, t, /) can be determined analytically. The result is summarized by the following 
theorem: 

Theorem 1 - IKT statistical model for CNSFE - Let us require that the IKT statistical model {/(x, i),r} 
satisfies Axioms #l-#3 and furthermore that: 

1) The CNSFE problem admits a smooth strong solution in T x I. 

2) The mean-field F(x,i;/) is defined as 



F(x, t- f) = FH(r, t) + F„(x, t; f) + F,(r, u, t; /) + u • VV+AF(x, t; /). 



(38) 



Here Fh is the Navier-Stokes acceleration, given by Eq. \142^ [see Appendix A] , F^ is the relative kinetic acceleration 
defined as 



-*/.v,.„„, "^+!^vin^^^ 



F„(x,i;/)-^Vlnp , 

2 2pi 



X' - 



with 



A 



D fpi 



D fpo + P 



+ nK{Y,t) 



■ Dt\p J '^ Dt \ p 

and with K{v,t) being prescribed by Ea. il 36]) in Appendix A. Finally, Fa(r, u, i;/) is defined as 

'Pi 



F,(r,u,t;/) = i[V-n-Vpi]+-^ 
p — 2pi 



V • Q - V In 



Q 



(39) 



(40) 



(41) 



3) F(x,i;/) is defined up to an arbitrary real gauge field AF(x, i;/) satisfying the gauge condition \21]) . We shall 
require that AF(x,i; /) is a smooth vector field analytic with respect to the Newtonian velocity vector v £ [/. 

4-) The BS entropy and the velocity moments J dvG{x,t)f{x,t) evaluated for {G'(x, i)} = {l, v, iy^, uu, iu^u} 
_ u 

exist for all (r, t) G fi x /. 

5) In fl X I the KDF /(x, t) admits the correspondence principle defined by Eqs.iSl]). 

6) Let us introduce the decomposition 



where 



P{f{t)) ^ I dv 
n 



^^S{f{t))^P,{f{t)) + ^^ST{t), 

Pi(/(i))^P(/(i))-|5T(i), 



'jp^A + pir,t)V.Y+'jp^{v.Q-VlJEl 
2pi 2pi [ \ p 



Q 



(42) 
(43) 

(44) 



and MrSTit) denotes the global thermodynamic entropy production rate [defined by (|162p ]. Then, we require that for 
all t £ I the pseudo-pressure Po{t) is determined so that 

Pi(/(i)) = 0. (45) 

It follows that: 

Til j The local Gaussian distribution function |J7| ) is a particular solution of the IKE H^ if and only if the fluid 
fields {Z} satisfy the CNSFE problem [see Appendix A]. For a generic KDF /(x, t), introducing the representation 



/(x,i) = /M(x,i)/i(x,i), 



(46) 



it follows that the reduced KDF h{x,t) satisfies the integral IKE 
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/i(x, i) = /iom,t„x) exp J - y di'^-|— • F,(x(t'), i'; /) i , 



(47) 



ho{'X.o) being a suitable initial KDF. 

TI2J In the case of a general non-Gaussian KDF /(x, t), the velocity-moment equations obtained by taking the 
weighted velocity integrals of Eg. jT^ with the weights {G(x,i)} = {l,v, u^/3} deliver identically the fluid equations 

TI3J For all t ^ I the pseudo-pressure po{t) must satisfy the ODE: 



d 



:Po 



dt'^^ I f dr^ 

2 J pi 



Sp{t)+Q{t) + -STit) 



(48) 



where 



Spit) 



2 J pi 

n 



dt 



(p - (j) + uT) + Y ■ \7 {p - (j) + nT) 



drpV ■ V, 



(49) 



Q{t) 



dv 



3^ 

2pi 



V • Q - V In 



Pi 



Q 



(50) 



TI4 j The BS entropy S (fit)) satisfies for all t £ I the H-theorem i35\) . Instead, for an isothermal fluid the constant 
H-theorem i36\} holds. 

TI5J In validity of Eq. (JS^ , for all (r, i) £ CI x I the kinetic pressure pi(r,t) is strictly positive. Furthermore, the 
initial value Po(^o) *s uniquely determined by prescribing the condition of vanishing of the initial BS entropy in the 
case f(to) = fM{to) [see Axiom #2, entropic principle]. 

Proof - First, it is immediate to prove that, in validity of Eqs. (|38|) - (|4T1) . f]yi{x,t) is a particular solution of the 
inverse kinetic equation (fT4|) . The proof is analogous to that given in Refs.[2, 5] and it follows by direct substitution 
of the distribution /M(x,i) in the same equation (Proposition Tli). This implies that /M(x,i) satisfies necessarily 
the integral Liouville equation ([T5|) . so that 



/Af (x, t) = fM{Tt,t,^, to) exp 



dt' 



d 



9v(t') 



F„(x(t'),t';/) 



d 



dr{t') 



V{r{t'),t') 



Therefore, in case of a non-Gaussian KDF the same equation manifestly implies also Ea. (|47|) . 

Instead, if we assume that in F x /, /(x, i) is a particular solution of the inverse kinetic equation, it follows that 
the fluid fields {p,V,p,T} are necessarily solutions of the CNSFE equations. This can be proved either in the case 
/ = /M(x,i) by direct substitution in Eq. p^ or, in the general case in which / 7^ fM{x,t) is an arbitrary smooth 
and strictly positive particular solution, by direct calculation of the velocity moments of the same equation, evaluated 
with respect to the weight-functions G{x,t) — (l,v, u^/3) (Proposition TI2). In fact, in validity of Axiom #2 and 
Eqs.([38 | - (|40p . the moments equations corresponding to Eqs. (|3T|) yield respectively: 



dp{r, t) 
dt 



V-[p(r,t)V(r,i)]-0, 



(51) 



^p(r, t)V(r, i) + V • [p{r, i)V(r, t)V(r, t)] = P;^V(r, t) = pFn, 



(52) 



d_ 
di 



Pi(r,t) + V-[V(r,i)pi(r,t) + Q] 



= A 



V-Q-Vln( ^ 1 -Q 



Q-Vln(^). 



(53) 
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The first two equations coincide, respectively, with the continuity and Navier-Stokes equations while the third one, 
thanks to Eq. (pn| . recovers the Fourier equation [see respectively Eqs. (|128p . (|129p and (|130p in Appendix A]. Let us 
now evaluate the entropy production rate -§^3 {f{t)) (Proposition TI3). First we notice that thanks to the Brillouin 
Lemma [H: 



which implies 



5(/(i)) = -yrfx/(x,t)ln/(x,i)< 
r 

< - /'dx/(x,i)ln/M(x,t) - - /'dx/M(x,i)ln/M(x,i) = 5(/mW), 
r r 

§-^Sifit)) EE J dx/(x, t)-^ ■ F(x, t; f) < 
r 

< g-^SifMit))^ J d^fMi^,t)—-Fi^,t-jM), 



(54) 



where respectively: 



^•F(x,i;/) = — yl+ — 
9v "■''■" 2pi 2pi 



d_ 
|;-F(x,,/m) = ^A + V.V. 



V • Q - V In ( ^ ) ■ Q 



v-v, 



(55) 

(56) 

(57) 



It follows that P{f(t)) < P{fM{t)), with P{f{t)) being defined by Eq. (|44l) . implying in turn that the inequality 

Q{t) > (58) 

is necessarily fulfilled for all f e /. Therefore, introducing the decomposition (|42|) - p3|) . the constraint ((45|) manifestly 
implies that Proposition TI3 must hold. In addition, the entropy production rate fulfills identically the constraint 



|5(/(,)) . p^ty 



(59) 



Hence, thanks to the entropy law (|14ip . necessarily the BS entropy S {f{t)) satisfies the H-theorem (|55|) (Proposition 
TI4). Let us now prove Proposition TI5, namely that the constraint (P5|) requires the kinetic pressure pi{r,t) to be 
strictly positive in i7 x /. For this purpose, we first consider the case / = /m and impose \/t € Iq = {t : t > to, Vt G /} 
that the constant-entropy condition S{fM{t)) — S{fM{to)) holds, requiring: 



|s(/„(,))./*A^^o. 



Then, thanks to the identity 
2pi 



dt 2j pi 2 J pi 

n n n 



— (p - (/) + nT) +V ■ \7 {p - (j) + nT) 



it follows that po (t) must satisfy the ODE 



dpojt) 
dt 



2 J PI 



(60) 



(61) 



(62) 



with Sp{t) being given by Ea.(|49|. Hence pi{r,t) is necessarily strictly positive in il x /. The same conclusion 
manifestly follows imposing instead 



d_ 

di 



Sifuit)) 



j'^lk^-ll'^^'^--'^ 



(63) 
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which imphes 



dpojt) 
dt 



fdr^ 
J pi 



Spit) + ^^Srit) 



(64) 



Analogous conclusion holds, thanks to the inequality ([58)) . also in the case / ^ /a/ [see Eq. pS)) ]. Finally, let us impose 
the initial condition for the initial kinetic pressure Po{to)- Denoting M = J drp, in view of Axiom #2 this requires 

n 



SifAiito))^^ I rfrplnpi+M 



In (27r) 



3/2 



drpln p = 0, 



which uniquely determines po (to). Q.E.D. 
Here it is worth noting that: 

• Eas. ([551) - (PT|) yield a realization of the mean-field F(x, t; /), of the type indicated above (see subsection 1.4) and 
holding for compressible NS thermofluids, which satisfies the CNSFE problem [see Appendix A]. 

• The expression of the vector field F(x, t; f) in Eas. ([551) - (|1T|) is determined, up to the gauge field AF, according 
to the form of the KDF as follows. In the case / — Jm it is obtained by solving explicitly for F the equation 



d d\ u 

+ v— ln/M-2— .F = 
at dr I vth 



av 



(65) 



and imposing the validity of CNSFE. The procedure is analogous to that outlined, for example, in Ref.[2|- 
Instead, in the general case in which / ^ /m, with / being a strictly positive KDF satisfying Axioms #l-#3, 
F(x, i;/) is determined by requiring it is of the form (|38I) - ([M)) with Fq(x, t;/) to be suitably prescribed. In 
particular, thanks to Axiom #3, Fa(x, t;/) is taken to be a polynomial of first degree in the relative velocity 

u. Therefore, it is necessarily of the form Fa(r, u,t;/) — Fi '{r,t;f) + uFa {r,t;f), with Fi '{r,t;f) and 
Fa (r,i;/) being respectively two suitable moments of the KDF /. Their precise form is obtained by imposing 
Axiom #1 and requiring that the velocity moments of IKE corresponding to the weight functions G{x,t) = 
(l, V, u^/3) coincide with CNSFE. For example, Fi '(r, t; f) follows by constructing the moment equation with 

respect to G = v. From IKE, utilizing the definitions for F/f(r,t) and F/f(x, f;/), it follows that Fi '{r,t;f) 
must be prescribed so that the equation 



p§-y + (V-n - Vpi) -pFh- J dv/(x, t)F(o) = 



(66) 



coincides with the NS equation. This yields the unique solution Fi \r,t;f) = - [V-II — Vpi] . Similarly, the 
expression for Fa (r, i; /) follows from the moment equation with respect to G = m^/3. This yields the solution 
Fi^'^ {r,t;f) = ^ V • Q - V In ( £i j • Q . The resulting expression for Fa(r, u, t; f) coincides with Eg.dlTt. 

• If / = /m(x, t) the extended fluid fields Q(r,i) and n(r,i) [see Eqs.(|32|)] vanish identically. As a consequence, 
in this case Fa(r, t) = 0. 



• Eqs. ([38l) -(|40 |) apply also in the case of an incompressible thermofluids [INSFE problem], and in particular for 
isothermal fluids [INSE problem; see Appendix A]. In the case / = /m(x, i) the corresponding mean- field 
F{x,t; Jm) is consistent with Refs.[2, |5|. However, in the case of non-Gaussian KDFs, the linearity condition 
here imposed as a kinetic closure condition on F{x,t;f) (see Axiom #3b) actually leads to a representation 
of the mean- field in terms of the extended fluid fields Q(r,i) and n(r,t) which is different from that adopted 
previously in Refs.@, Q- 

In addition, in agreement with the GENERIC dynamical model [lllli^lj it is possible to show that: 

• Both in the case of CNSFE and INSFE the IKT-statistical description permits to represent the reduced set of 
fluid fields {Zi} [see Eq.([7])] in terms of velocity and phase-space moments (bundle structure on T). 
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• The functional setting for /(x, i) can be suitably prescribed so that {/, F} is compatible with the physical 
observables [llj. In other words, the predictions for the fluid fields obtained in this way in terms of the statistical 
model are in agreement with the experimental observations. This includes, in particular, the prescription of 
appropriate kinetic boundary and initial conditions [5| . 

• The NS-DS is a deterministic, non-conservative, irreversible and non-canonical dynamical system. In fact, 
through the fluid fields, the mean- field defined by Ea. (l38p becomes generally explicitly time-dependent, while 
generally ^ • F 7^ 0. Hence, the NS-DS [see Eq.([5])] defined in terms of F is manifestly irreversible and 
non-conservative. Furthermore, it is possible to show that Eq. psp is non- variational, and hence intrinsically 
non-symplectic, so that it cannot be cast in local canonical form. 

V. DYNAMICS OF NS IDEAL TRACER PARTICLES IN DETERMINISTIC FLUIDS 

A basic consequence of the previous theorem is that the functional form of AF remains non-unique for ITPs. Its 
possible unique prescription requires, therefore, the adoption of suitable additional kinetic closure conditions. For 
this purpose in this section we intend to show that for all TTPs AF can be uniquely determined in such a way to 
fulfill the requirement of the Gedanken experiment (GDE-requirements #l-#6). 

A. Preliminary Lemma 

Let us first show that due to the kinetic constraints (|27p and (I^Hl) for TTPs the pseudo- vector i},{r,t) entering the 
evolution equation for n(r, i) [see Eg. ([^5)) ] is actually uniquely determined. In fact the following result holds. 

LEMMA to THM.2 - General form of n{r,t) 

If fli = fl X I denotes the existence domain of the fluid fields {Z}, let us assume that in the subset Qi C f2i in 
which |Vpi| ^ 0.' 

1) the real unit vectors n(r, i), b(r, i) and the pseudo-vector n{r,t) are all differentiate and suitably smooth; 

2) the unit vector n(r, t) satisfies the kinetic constraint \2T^ and the initial-value problem V2^] : 

3) moreover, Jl(r, i) satisfies the constraint I129\) and is defined in the limit pi{r,t) — > 0^ and also for arbitrary 
finite values of the kinetic pressure pi(r,t). 

It follows that, in the domain fli, Jl(r, i) necessarily takes the form: 



n(r,t) = b(r,t)x^^M + c(r,t)b(r,i), (67) 



where the GDE-requirement #6 implies 



c(r. i) = b X n • -^ = ftCr, t) ■ b(r, t) = -^(r, t) ■ b(r, i), (68) 

dt 



dn{T,t) 



with ^(r, t) denoting the local fluid vorticity. 

Proof - In fact, by definition |n(r, i)| = 1 and hence ''"^^'''' • n(r,t) = 0, so that there must exist a pseudo-vector 
Jl(r, i) such that Eq. ipS)) holds identically for arbitrary initial condition n{r{to),to) = n(ro,io)- Let us now impose 
the validity of the kinetic constraint (|27)) , implying 

-n(r, t) ■ b(r, t) = -n(r, t) ■ -b(r, t). (69) 

Hence, necessarily r2(r,i) satisfies the equation 

n(r, t) ■ ^b(r, t) = - [fl(r, t) x n(r, t)] ■ b(r, i), (70) 

i.e., due to the arbitrariness of the unit vector n(r,t), h{r,t) x n{r,t) = — ^b(r,f). This yields for 0(r, i) a general 
solution of the form 

f2(r, t) = b(r, t) X ^^M + c(r, t)h{r, t). (71) 

dt 

Substituting this solution in Ea. (P5|) and taking the scalar product of the resulting equation by b x n, in validity of 
GDE-requirement #6 Eq. (l551) follows identically, which uniquely determines the form of r2(r, i). 
Q.E.D. 
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Construction of TTP solutions 



Let us now prove the existence of the TTPs, particular solutions of the initial- value problem ([4]). More precisely, we 
intend to prove that for all t e / and an appropriate choice of the mean-field F + AF, the NS-DS (I2|) [or equivalent 
RD-NS-DS) ^], necessarily maps an arbitrary TTP initial state x(to) = Xq into a TTP state x(i) — Tt^ tX,,. For this 
purpose we impose that, consistent with GDE-requirements #l-#6, F is defined by Eqs. (|39|) - (|40|) . We intend to show 
that, as a consequence, for all TTPs both F and the gauge field AF [see Eq. ipT]) ] are necessarily uniquely determined. 
On the other hand, in view of GDE, one expects that the vector field F which characterizes TTP dynamics should 
not depend on the form of the KDF. If true, the result would clearly be conceptually important because it would 
imply the uniqueness of TTP dynamics in all cases. For reference, let us first consider the case of the Gaussian KDF 
/ = Im, leaving the extension to a non-Gaussian KDF to the discussion below. Then, the following theorem holds. 

THEOREM 2 - Existence and uniqueness of TTP dynamics 

In validity of THM.l, let us require that AF and F^ are analytic functions in T which are defined also in the limit 
Pi(r,i) — > 0"*". Then it follows that: 

T2i) The initial-value problem defined by (01) admits particular solutions fulfilling the GDE-requirements #l-#6, 
here denoted as TTPs. In particular, for all t G I they are characterized by a relative kinetic velocity defined by 
Eq. Ii23\) and is such that: 

A) the local magnitude of the relative kinetic velocity |u(t)| is determined by the equation: 



u(i)| = uth{r,t) = I3vth{r,t), 



(72) 



with /3 > being independent of (r,i); 

B) n(r, t) satisfies both the constraint equation 

C) ft{r,t) is uniquely determined by Eqs. |g7p and ^ 
T22) For arbitrary TTPs, the mean-field F{r, Uth, t; fM 



d the initial-value problem 



AF has necessarily the unique representation 



F(r, utk,t; fM) +AF = Fh + Uth ■ VV 



^-ln(pi; 
2 Dt ^^ ' 



- I3uthftir,t) X n. 



(73) 



This implies that, in the case of the Gausssian KDF f 
AF(r,Uth,t) has the unique representation 



fM, o,nd for arbitrary (3 > 0, the vector field AF 



AF = AFo + AFi, 

1^ V In p + 4^ V In (pi) (/32 - i) 



AFo{r,Uth,t) ^ 

AFi(r, uth,t) = I3uth^{r, t) x n. 

T23) Particular solutions of the form i23\) which fulfill requirements A-C must satisfy the initial conditions 



(74) 



r(to) = To, 
u(io) = /3vth{ro,to)n{ro,to), 



(75) 



with /3 > being an arbitrary real constant, pi(ro,to) the initial kinetic pressure and n(ro,to) a unit vector satisfying 
the orthogonality condition 



n{r o, to) • Vpi(ro,to) =0. 



(76) 



Proof - T2i— T22) For generality let us assume that |Vpi| ^ everywhere in fl x I. Then, it is sufficient to prove 
the theorem in the subset f2i C fli = fl x I in which |Vpi| ^ 0. Let us show that in ili, for an arbitrary non-negative 
constant (3 S M"*" , a particular solution of the initial- value problem Q of the type p3|) , which satisfies requirements 
A-C, exists and is unique. In fact, let us assume that |u(i)| is of the form (I72|) . with /3 = /3{r, t) > denoting now an 
arbitrary smooth real function of (r, t) defined in fi x J. It is immediate to show that necessarily 13 must be everywhere 
constant with respect to (r,i) infl x I. Indeed, Eq.(|3]) requires 



^ i(3vth{r,t)) =. ^n-Vlnp+^Pvthir,t)^lnip,)- 
+ ^n ■ Vln (pi) (/32 - 1) -)- n ■ AF, 



TM^^t) = ^ [^Vlnp+^-fVln{p,) (/?2 _ r 
+ ^AF-{l-nn)^nir,t)xn 



1 nn) 



(77) 
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On the other hand, imposing the constraint ((27)) requires necessarily, thanks to the Lemma, that n(r, t) must satisfy 
the initial- value problem (P5|) . We require for consistency 

n-Vln(pi) = 0, (78) 

v'^ — 

^n- Vlnp + n- AF = 0. (79) 

Hence, since by assumption AF is defined also in the limit pi{r, t) — > 0+ (or equivalently f3 — )► 0+), it follows necessarily 
that 

AF = AFo(r, uth,t) + Ari(r, uth,t), (80) 

AFi(r,Ut,„i)-(i-nn) =/3uthn(r,i) xn(r,t), (81) 

with J7(r,t) being given by the Lemma and 

d/3(r,i) 



dt 



= 0. (82) 



Therefore, due to the arbitrariness of v = V(r, t) + uth{j^, t)n{r, t), (3 is necessarily independent of (r, i); furthermore, 
AF(r, Ut/i, t; /m) and F{r,Uth,t; fm) are necessarily of the form (|74p and ([73l) . T23) Finally, the initial conditions 
(|75p are an immediate consequence of Ea. (|23l) and the requirements A-C. 
Q.E.D. 

Let us now consider the extension of the theorem to the case of a non-Gaussian KDF. We notice that the gauge 
field AF evaluated for the state of a generic TTP, can always be identified with the vector field 

AF = AFo + AFi-Fa(r,Utft,i;/), (83) 

where AFq and AFi are still given by Eq. ([74|) . while Fa(r, Uth,t; f) is prescribed according to Eq. (j4T|) and computed 
for u = Uf/j. As a consequence, it is immediate to show that, for TTPs, this prescription of AF warrants the uniqueness 
of the vector field F(r, Uth, t; /m), in agreement with GDE. Therefore, for TTPs, its form is independent of the form 
of the KDF, namely F(r,Ut,i,i; /m) =F{r,Uth,t; f). 

C. Implications and physical interpretation 

Let us briefly analyze the implications of THM.2. First, we remark that by construction for all TTPs the mean-field 
acceleration F -I- AF is unique and independent of the form of the KDF f{x,t). As a consequence TTP particular 
solutions [of the initial- value problem Q] realize a classical dynamical system with existence domain Fi x J, Fi denoting 
a suitable subset of the phase-space F (see related discussion in Section 6). This is defined by a homeomorphism of 

the form T^^JJ^^ : yo ^ y{t) = t'^^J^^Yo, where for all t £ /, 

y(i) = {r,u(i)=u,„(t)} (84) 

and uth{t) is given by Eq.([23l). Manifestly t'^'^J^^ is a subset of the RD-NS-DS defined by Eq.® [or equivalent of 
the NS-DS defined by Eq.©]. 

Furthermore, let us assume that the initial conditions for Eq.Q are of the form y(to) = {yq, u^ = VLth{to)}: where 
Ut/i(io) is prescribed by Eg. d^^ while the initial unit vector n(r,io) satisfies the constraint (PS)) at r = Tq. Then, 
thanks to THM.2, it follows that for all t € I, y{t) it is necessarily of the form (IMl) . i.e., it defines, for all t E I, a. 
TTP. In addition, thanks to Eg . ([5^ it follows that /3 is constant and therefore is uniquely determined for each TTP 
by the initial state y{to), namely /? = - — 1"° '^ y 

We notice that, by construction, TTPs are uniquely associated to the local state of the fluid. As a consequence, this 
permits to determine also for the remaining NS-ITPs an explicit representation of the mean- field which is necessarily of 
the form F(r, u, t; /), namely it depends explicitly on the KDF. Similarly the gauge-field is of the type AF(r, u, t; /). 
Both hold for arbitrary u G [/ and satisfy at the same time the requirements posed by THM.2, namely that for 
u = Uth, the sum of the two vectors F(r, u, t; f) + AF(r, u, t; /) must reduce to Eg. ([75)1 . For definiteness, in validity 
of THM.l, let us consider the case of a generally non-Gaussian KDF /(r, u, i), with /(r, u, i) denoting a smooth, 
strictly positive function which satisfies at the same time Axioms #l-#6, i.e., is a particular solution of Eq. (fT4l) . Let 
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us therefore determine AF(r,u,i;/) in such a way to fulfiU the constraint equations ((73|) and the gauge condition 
(PTj) . For this purpose it is sufBcicnt to let 

AF{r,u,t;f)^AF{r,nt,„t)lp^^, (85) 

with AF{r,Uth,t) being defined by Eg.dTil). 

Finally, an interesting issue concerns the physical interpretation of the evolution equation for the unit vector n(r. t) 
[i.e., the direction of the particle relative velocity] and related pseudo-vector il. In fact, that Eas. ((28l) are similar 
to the Euler equations for a rigid body rotating with angular velocity il^ = —fl. This suggests that due to the 
GDE-requirements #5 and #6, two different physical effects contribute to il. These are due both to the rotation of 
the unit vector n(r,i) as determined by Ea. ((28| as a consequence of non- uniform specific kinetic pressure, and the 
contribution of fluid vorticity specified by Eg. d^ . Indeed, from Eq. ([7T|) . denoting by — ^^ = ^ + V(r,i) • V the 
fluid convective derivative and since by construction Vu — — VV, it follows 

dh(r f) Dh(r f) 
nr= ^= \' ' X b(r, t) - c(r, i)b(r, t) = ^-^ x b(r, i) + (u • V) b(r, t) x b(r, t) - c(r, i)b(r, t), (86) 



dt ' ' ' V ' / ^ w ^^ 



where 



(u • V) b(r, t) X b(r, t) - c(r, t)b(r, t) = (87) 

= e - 7^ [b X V(Vpi • V) + b X (Vpi • V) V] . 

In the last equation the first term on the r.h.s. denotes the vorticity ^ = V x V(r, t). This means that near a vortex 
the motion of TTPs is qualitatively similar to that of a rotating rigid body. However, by inspection of the remaining 
terms in Eas. ([5Sl) and (|57|. it is evident that more complex particle-acceleration effects may be present, which are 
driven by time-dependent pressure and velocity-gradients contributions. 

VI. THE TTP-STATISTICAL DESCRIPTION 

In this section we introduce a statistical description associated to the ensemble of TTPs, denoted as TTP-statistical 
model, which is represented by the couple {/i,ri}. We intend to show that, like the IKT-statistical model {/, F}, 
also {/i,Fi} determines uniquely the time evolution of the complete set of fiuid fields {Z}. However, the result is 
conceptually important because Fi is a reduced- dimension subset of T . For this purpose, we notice that if xith denotes 
the relative velocity of an arbitrary TTP endowed with a relative-Newtonian state y = (r, Ut/;), then \ith spans the 
subset of velocity space U: 

C/i = { u| u e [/, u = uth = Pvth{Y, t)n{r, t), n(r, t) ■ b(r, t)^0,(3eR+}. (88) 

2 

Here by construction /3^ = ^""' , is a constant independent of (r,i). As a consequence, it follows that Fi is the 

subset of F, Fl = ri X Ui, with UiCU = R^ and dim{Ui) = 2. 

To define a KDF /i(r, u^^j, t) on Fi let us first consider the PDF defined on U in terms of the KDF /(r, u, t). For 
a prescribed IKT-statistical model {/, F}, the corresponding conditional velocity PDF on Ui is defined as 

^^^'^'"^^'^^= /rfry/(r,u,0 - ^^^^ 

Here, introducing for Uj^ a representation in terms of the spherical coordinates {u = Uth,'^T'&) and requiring the 
Uz = b(r,i) it follows by construction that both the conditional PDF and KDF defined by Eqs. (j89|) and (|90|) are 
independent of the angle i?. In particular it follows that by definition J drif{r, u,t) = J dvS {-d — 7r/2) /(r, u, t) and 

Ui u 

hence dr/ = u^dudip, while the corresponding phase-space measure is dxi = d'^rdr/. Thus, the conditional KDF on Fi 
is defined as 

h{r,uth,t)^p{r,t)h{r,uth,t). (90) 



In particular, if f{r,Utfi,t) coincides with the Gaussian KDF ((TT)) . it follows that 

2pir,t) 



fi{r,uth,t) 



7r3/2^;3^ 



exp{-/?^} = /iM(r, ut^,i). 
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(91) 



The main result can be summarized by the following theorem. 

THEOREM 3 - TTP-statistical model for CNSFE 

Let us require that the IKT- statistical model satisfies THMs. 1 and 2. Then it follows that the conditional KDF 
fi(t) = fi{r,Uth,t) defined by Eas. \89]) and Ii90\) has the following properties: 

T3i) It is a particular solution of IKE which is independent of the angle i? 
emph. 

T32) It satisfies the functional constraints: 



Sifiit)) 



J dr]fi{v,uth,t) ^ p{T,t), 

Ui 

jrh) I dv^h (r, Uf h , i) = V(r, i) , 
/ dri\u'^fi{r,Uth,t) =pi(r,t), 



(92) 



/ dxi /i (r , VLth , t) In /i (r , Uth ,t) ^ St {t) 



(correspondence principle) . 

TSa) Its velocity moment equations, determined from IKE in terms of the weight-functions {G(x, t)} — {l, v, gW^}, 
imply again Eqs. i31\) and therefore, together with Ii89\) and the constraint equation ^8^ , they coincide again with 
CNSFE, so that in particular /i (i) satisfies the entropy law. 

Proof - T3i) In fact, due to the hypothesis and the definitions ([M)) and (lUU]) . it follows that in Ui 

fi{r,uth,t)=2fir,uth,t), (93) 

and is therefore by construction independent of the angle ■&. Hence, in the subset Fi x /, fi{r,Uth,t) is a particular 
solution of IKE, in the sense that it satisfies by construction the statistical equation 



d d 

— /i (r, ut/i, i) + V • V/i (r, Uth , i) + F(r , Ut^ , t) • /i (r, Uth , t) 

at duth 



(94) 



+fi{r,uth,t) 



d_ 
dv 



F(r,u,i,/)/(r,u,t) 



= 0. 



T32) The proof follows by noting that, thanks to Eq. (j93|) and the fact that f{r,Uth,t) satisfies by construction 
the correspondence principle (see THM.l), the conditional KDF fi{r,Utfi,t) fulfills identically the correspondence 
principle (|92|) too. 

T33) Due to T3i),T32) and Eg. dMl) the moment equations coincide necessarily with CNSFE and hence, in particular, 
consistent with the second principle of thermodynamics [see Eq. (|13ip ]. fi{t) satisfies the weak H-theorem 



|5(/iW)>o. 



(95) 



Q.E.D. 

Wc remark that here: 

• /i(r, Uth, t) denotes generally a non- Gaussian conditional KDF. A particular solution is provided by the Gaussian 
conditional KDF defined by Eq. ipTj) . 

• /i(r, Ut;i,t) satisfies the inverse kinetic equation ([M]) . 

• {/i,ri} is a reduced-dimension statistical model for CNSFE problem. The fluid fields {Zi} are uniquely 
advanced in time by means of Ea. (|94p . or equivalent, by means of the integral Lagrangian IKE (jlSp . 

• fi{r,Uth,t) determines uniquely the time-evolution of the fluid flelds {Zi{r,t)}. 

• In view of the discussion presented above after THM.l, {/i,ri} applies also to incompressible fluids described 
either by the INSFE or INSE problems. 



In the following we analyze basic implications of this result. 
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VII. TTP-DYNAMICS IN STOCHASTIC FLUIDS 

A remarkable application of the TTP dynamics concerns the modelling of stochastic tracer-particle dynamics in 
stochastic fluids, such as for example due to temperature and pressure fluctuations (thermal fluctuations). Thermal 
fluctuations are important in a wide variety of mesoscopic flows (see for example Refs. |29l - [3l| ). Theoretically, they are 
usually treated within the framework of fluctuating hydrodynamics, an approach pioneered by Landau and Lifshitz 
[32,[33- In this framework extra (stochastic) terms are added to the fluid equations to model possible stochastic effects, 
so that generally the functional form of the corresponding fluid equations is actually modified with respect to the 
customary fluid equations. For example, in this case the latter may typically include higher-order spatial derivatives of 
the fluid fields. The numerical solution of fluctuating hydrodynamic equations may present, as a consequence, serious 
difficulties (which are nevertheless also present in the case of the incompressible NS equations for isothermal fluids) . 
A possible alternative is represented by particle simulation methods based on kinetic theory. In such an approach, 
unlike fluctuating hydrodynamics: 

1) Requirement ^1: the functional form of the corresponding fluid equations and of the related initial-boundary 
value problem is left unchanged |9l. llCt l34|. i.e., the differential operators appearing in the fluid stochastic equations 
are the same ones entering the customary fluid equations in the absence of stochasticity. 

2) Requirement #2: the stochastic fluid fields are assumed to be strong solutions of the stochastic CNSFE problem 
(see Appendix A). 

An approach of this type can be achieved by means of the TTP-statistical model. A convenient representation of 
stochastic fluid fields of this type, and fulfilling Requirements #1 and #2, is provided by Eg. ((30|) . In this case the 
fluid flelds are assumed to depend on a suitable set of stochastic variables a = {ai, i = l,k} G Va Q R'^, with fc > 1, 
by assumption all independent of (r,t) and endowed with a stochastic probability density g{r, t, a) on Va [^jllOilSJI- -"-^ 
must be remarked that "a priori" the parameters a = (ai, ..., a„) and the related probability density g{r,t, a) can be 
set arbitrarily. Thus, they can in principle be chosen to provide prescribed mathematical models of stochasticity. In 
the following we shall assume in particular that the a's are also independent of (r,t). As a consequence, introducing 
the stochastic averaging operator (|172p , the fluid fields can be represented in terms of the stochastic decomposition 

Z(r, t, a) = {Z{r, t))^ + SZ{v, t, a), (96) 

{Z{r,t))^ and SZ{r,t,a) denoting respectively the corresponding stochastic-averages and stochastic fiuctuations of 
the fiuid flelds. 

In particular, in contrast to fluctuating hydrodynamics, in the present approach the functional form of the fluid 
equations is left unchanged. It follows that the precise form of the stochastic-averaged fluid fields {Z{r,t))^ and of 
the stochastic fluctuations SZ{r,t,a) depends solely on the model of stochasticity adopted, i.e., the choice of the set 
{(?(r, t, a), Va}- This means that its realization may generally depend on the possible sources of stochasticity adopted, 
namely: 1) Stochastic initial conditions: in this case the initial fluid fields Z(r,to) = Zo{r) are assumed stochastic, i.e. 
of the form, Zo{r, a) — {Zo{r, a))^+SZo{v, a), with {Zo{v, a))„ and dZo{r, a) being suitable vector fields. 2) Stochastic 
boundary conditions: this occurs if the boundary fiuid fields Z^ (r , f ) | g^ are prescribed in terms of a suitable stochastic 
vector field of the form Zt^:{r,t,a)a\sQ = {Zwi'r,t,ct))^\g^ + i5Zu,(r,i, a)|^Q. 3) Stochastic forcing: in this case the 
volume force density acting on the fiuid is assumed stochastic, i.e., of the form f(r, i, a) = (f(r, i, a))^ + Si{r,t,a), 
being (f(r,i, a))^ and Si{r,t,a) suitable vector fields. 

A. Langevin dynamics in fluctuating fluids 

As stated above, a fundamental consequence of TIIM.2 is the uniqueness of the deterministic equations of motion 
for arbitrary TTPs belonging to a compressible or incompressible, thermal or isothermal NS fluid. It is immediate 
to show that the dynamics of TTPs is unique also when the same fluids are considered stochastic. It follows that 
the relative state y of a generic TTP advances in time by means of a stochastic dynamical system (DS), namely 
the flow generated by the initial value problem associated to the stochastic equations of motion for TTPs (Langevin 
equations): 

^ = Pvth(r, t, a)n(r, t, a) + V(r, i, a), 

^=F„(r,/3«,^n(r,t),i,a), (97) 

y(io,a) =yo(a), 

where the unit vector n(r,i,a) satisfies the stochastic initial- value problem 

''"'dt*^°^ = ^(r, t, a) X n(r, i, a), ,^^. 

n(ro,to,a) = no(r,a). 



20 

In particular, it follows that the stochastic mean-field F is provided by Eq. ((38)) . with F^ = Fu{r, /3vthn{r,t),t,a) 
being identified with 

F„(r, l3vt,Mr, t),t, a) - l3vth{r, t, a)n{r, t, a) x n(r, t, a) + /3n(r, t, ^^E^^lh^^L^^ (99) 



where n(r,i,a) is given by Eq. ipS)) (see also Lemma to THM.2). The flow associated to the Eqs. d^Tl) and 

Tt^A:yo^yit.a)^Tt^,tyo, (100) 

is referred to as stochastic TTP dynamical system (TTP-DS). Here Tt^^t is a measure-preserving evolution operator 
associated to the Newtonian vector field X(x, i, a) = {v, F(x, i, a)} and F(x, t, a) the stochastic mean-field satisfying 
the initial- value problem (l99l) . Thus, by definition, the TTP-DS is uniquely prescribed by the instantaneous state of 
a generic TTP x(i,a) = x = (r, v). Therefore, we conclude that: 

• The initial-value problem ([W)) can be viewed as a stochastic model of particle motion in compress- 
ible/incompressible NS fluids, describing the dynamics of TTPs in stochastic thermofluids. 

• Based on the TTP-statistical model developed in the previous section (see THM.3) the time-evolution of the 
stochastic KDF /i(r, Uf/j, i, a) is uniquely prescribed. In particular, independent of the choice of the stochastic 
model {g(r, t, a), Vq}, the dynamics prescribed by (jlOOp preserves the exact form of the stochastic fluid equations 
and is unique. 

• The stochastic KDF /i(r, Ut;i, t, a) determined by Ea. ((T5|) is necessarily a particular solution of Eq. ([M)) . 

• The stochastic KDF /i(r, Ut;j,i, a) prescribes uniquely the time-evolution of the stochastic fluid flelds 
{Z(r,t,a)}. 

B. Fokker-Planck description in strong turbulence 

Let now analyze the time evolution of the stochastic KDF /i(r, Ut/j,i,a) represented in terms of the stochastic 
decomposition 

/i(r,ut,,,t,a) = {}i{Y,nth,t,a))^ + Sfi{r,Uth,t,a), (101) 

with {fi{r,Uth,t,a))^ denoting the stochastic-average defined by Ea. (jl72p . Then, requiring that the stochastic PDF 
g{r,t,a) is homogeneous and stationary [see Appendix B], i.e., that g = 5(a), it is immediate to obtain from IKE 
[see Eq. p^ ] the stochastic statistical equations advancing in time {fi{r,Uth,t,a))^ and Sfi{r,Uth,t,a). These are 
explicitly 

(i)a (/l). = - iSLSh)^ EE (C)„ , (102) 

(L)^ Sh = -5L {(/i)„ + Sh} + {SLSh)^ , (103) 

where the streaming operator L has been similarly represented as i = (i)^ + 6L. Eqs. (ll02p and (|103p are formally 
similar to the Vlasov equation arising in the kinetic theory of quasi-linear and strong turbulence for Vlasov-Poisson 
plasmas |34l - l37| . As is well-known, the construction of the precise form of the operator (C) appearing in the stochastic- 
averaged kinetic equation [i.e., Ea. (jl02p ] represents a task of formidable difficulty. The reason is that it requires 
constructing a formal perturbative solution of the equation (I103P for the stochastic perturbation 6f. To obtain a 
convergent perturbative theory, however, this usually requires the adoption of a suitable renormalization scheme in 
order to obtain a consistent statistical (kinetic) equation for (/i)^, [see earlier approaches developed in Refs. |38l - l40| 
which pertain to the statistical treatment of particle dynamics and the specific application to Vlasov-Poisson plasmas] . 
Nonetheless, in the case of weak-turbulence, the stochastic-averaged kinetic equation [i.e., Eq. (|102p ] is known to be 
amenable to an approximate Fokker-Planck kinetic equation advancing in time (/i)„ alone. Analogous suggestions 
are provided by p henomenologically-based Markovian Fokker-Planck models of small-scale fluid turbulence [see for 
example Refs.dlHil]. 

This raises the issue of (the construction of) a possible representation of this type for the stochastic-averaged 
operator (C)^ which has the following properties: 

• Property ^1: it holds at least locally in the velocity space f7i, in a suitable space to be specified. 
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• Property ^2: it holds in the case of "strong turbulence^\ namely when there are fluctuating quantities such that 
their stochastic fluctuation is comparable in order of magnitude with their corresponding stochastic averages. In 
particular, denoting by C a dimensionless infinitesimal parameter, in the following the strong turbulence regime 
is defined in such a way that [56[ 

I^V|^|(V)J[1 + 0(C)], (104) 

\Spi\^\{p,)J[l + OiO]. (105) 

The two requirements are mutually consistent and, for arbitrary choices of the velocity stochastic fluctuations 
(5V, are required by the Navier-Stokes equation. 

• Property ^3: it is applicable also in the case in which (/i)^ is generally non- Gaussian. 

• Property ^4: it does not rely on renormalization theory. 

• Property #5: the KDF is assumed of the form /i = /i(r, C^^^Uth, t, a), namely it exhibits slow dependence with 
respect to the velocity Uth- 

Notice that the previous properties are assumed to hold for arbitrary smooth non-Gaussian KDFs. In the specific 
case of a Gaussian KDF Property #5 requires necessarily that -^^7-^ ^ 0(C^'^). For TTPs this implies ^j^ — 

/3n ^ O (C^^^)j namely (3 ^ O (C^^^) and also "*'' )j""''^° ^ O {(). Therefore, the required slow velocity dependence 

th 

effectively limits the validity of strong turbulence theory to the subset of velocity space in which such an ordering 
holds. 

Regarding, in particular, the form of the KDF /i(r, C^/^Uj^, t, a), it must be noted that, even if the latter coincides 
locally with a Gaussian KDF /7\/i(r, (^^/^Ut^,, t, a), its stochastic average (fi)^ still remains generally non-Gaussian. 
Here we intend to prove that an explicit representation of (C)^ fulfilling the previous properties can be determined, 
based on the IKT approach earlier pointed out in Ref.flOl. In this case, the previous problem is exactly solvable 
provided: 

1. /i depends on the stochastic variables a only through the fluid fields, i.e., p{r,t,a), pi{r,t,a) and V(r,i,a), 
and hence also uth = ^Vthi'^i t, a) and the unit vector n = n (r, t, a), namely is of the form 

/iEE/i(r,Ci/2ui,„t,p,pi). (106) 

2. /i is generally a non-Gaussian KDF which is analytic in Uth, p and pi. 

3. The kinetic pressure pi and the pseudo-pressure po = (po)^ satisfy the ordering 

Pl^P^^OiC). (107) 

Po 

4. The mass-density perturbations are weak, in the sense 

M 

5. The stochastic fluctuation of the KDF Sfi{r, (^'^Uth, t, a) is such that 

(5/i(r,Ci/2ut,j,i,a 



0(C). (108) 



(/i(r,Ci/2u,,,i,a))^ 



0(C). (109) 



Invoking for Suth the representation Suth = Sp g,"t '' + Spi ^^^ — 6Y, as a consequence of the previous requirements. 



/i can be Taylor-expanded with respect to Sp, Spi and Suth yielding 

/i(r,C"'u,i,i,p,pi) = A(r,((;"2u„,)^,i,(p)„.(pi)„)+ 



t^l'" 



n 

.11=1 



{'"m: ^ *■!: -^ ""'• - 1) ''^'■''"''■'- <">• ■ <">-' 



(110) 
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Due to the previous ordering assumptions and the requirement of analyticity, the series converges uniformly in Fi 
and permits the explicit determination of both (^fi{rX^^^Uth,t,pi)) and Sfi{r,(^^'^Uth,t,pi). As a result, after 
straightforward algebra the operator (C)^ becomes explicitly [lOl 



(C)„ = Cfp ee 



'.,ji',A;— 1 " 



(111) 



where C^.j^fc are the tensor Fokker-Planck (or Kramers-Moyal) coefficients 

c,,,^k^-^{sFiSpy{dpiyiduu,)') , (112) 



(E) denotes the tensor product and the summation is carried out on i,j,l from to oo, with n = i + j + k > 2. 
Therefore (C)^, takes the form of a generalized Fokker-Planck (F-P) operator (Cfp)- Remarkably, Eq. ljllip satisfies 
by construction Properties #2-#5. 

As a final point, we comment on the implications of the GDE-requirement #6 which concern TTP dynamics and 
hold in validity of the strong turbulence formulation developed here. In fact, the kinetic constraint (P^ implies that 
time evolution of the TTP relative velocity depends on the stochastic fluctuations of both kinetic pressure Spi and 
fluid velocity (5V. Qualitatively this means that, even in the case in which pressure fluctuations are negligible, TTP 
dynamics can still exhibit turbulent motion through the fluctuations S'V entering in the fluid vorticity. 

C. Statistical irreversibility 

We first notice that due to the integral IKE dTSl) . if /i(io) is strictly positive then, for all t > to m the existence time 
interval /, also fi{t) is necessarily so. This manifestly implies, in turn, that in such a case the stochastic-averaged 
KDF (/i(x, t))^ = {fi{r,Uth,t,a))^, solution of the stochastic-averaged statistical equation (I102p . is necessarily 
strictly positive. It is immediate to show that, thanks to THMs.l and 3, {fi{t))^ must satisfy in Fi x / also a weak 
H-theorem of the form: 

dsiiMM > 0. (113) 

dt - ^ ' 

For definiteness, let us assume that at the initial time to both /i(x,t) = fi{r,Uth,t,a) and (/i(x,t))^ are strictly 
positive and admit the BS entropies S'(/i(x,to)) and S'((/i(x,io))Q,)- Then, we notice that if g{x,t) is an arbitrary 
strictly positive function such that / d'x.ig{K,t) — fJ.{ft), the majorization 

y"rfxi/i(x,i)ln/i(x,i)>y dxi/i(x,i)lng(x,i) (114) 

necessarily holds (Brillouin Lemma 28]). On the other hand, assuming that the stochastic-averaging operator (•)^ 
commutes with the phase-space integral operator Jp dxi, 



it follows 



This yields in turn 



y dxi(/i(x,t))„ = /y dxi/i(x,i)\ , (115) 

/ dxi/i(x,i)ln/i(x,i)> / dxi/i(x,t)ln(/i(x,t))„. (116) 

[ dxiA(x,t)ln/i(x,t)\ > / dxi(A(x,t))„ln(A(x,t))„. 



(117) 

/ a -^Ti 

The last inequality implies manifestly that 

5((/i(x,i))J>(5(/i(x,i))„. (118) 

Therefore, in view of the entropy inequality (pS)) . (/i(t))^ satisfies necessarily the weak H-theorem (|113p . This assures 
that for all t > to in the time interval /, {fi{t)) admits the BS entropy integral S {{fi{t))), i.e., that S {{fi{t))) is 
defined for all t > to in I. As a consequence (/i(i)) exhibits an irreversible behavior. 
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D. Comparisons with the HRE statistical model and solution of the Closure Problem 

In this section we display the relationship between the statistical models obtained here - i.e., both the IKT and 
TTP statistical models {/, F} and {fi,Ti}, prescribed respectively by means of THMs 1 and 3 - and the customary 
statistical treatment of turbulence due to Hopf, Rosen and Edwards {HRE approach, [24| - |26j : see also |44| - |48j : for a 
review see [43, [50|)- Since the latter is usually developed in the case of incompressible fluids, we shall restrict the 
following analysis to such a case. 

The HRE approach, which in its original form applie s only to incompressible and isothermal NS fluids, is based on 
the introduction of a suitable statistical model [2J-|26[ {fn,^}, here referred to as the HRE statistical model, with 
T = n X U. In this case, denoting by u = u(t) = v — V(r, t), fn is identified with the velocity PDF 

/ff(r,u,i)=p„J(u) (119) 

{HRE velocity KDF), po denoting the constant mass density characterizing an incompressible fluid. Upon identifying 
the mean-field F(x,i) with the fluid acceleration F/f (r,i) [see definition given by Eq. (jl42|) in Appendix A], it follows 
that by construction fn is a particular solution of Eq. p^ . As a consequence, its velocity moment equations, evaluated 
with respect to the weight functions {G} = {l,v}, coincide respectively with the continuity and the NS equations 
(I128p . (|129|) (see Appendix A). It is important to remark, instead, that by construction the fluid pressure p{r, t) cannot 
be determined as velocity moment of fif{r,u,t) of the form J dvG{r,v,t)fif{r,u,t) in terms of a weight function 
G(r, V, t) which is independent of the same fluid field. Hence, in the HRE approach the fiuid pressure has to be 
suitably prescribed. Thus, for incompressible NS fiuids it is identified with a solution of the boundary- value problem 
associated to the corresponding Poisson equation. 

Let us now pose the problem of the connection existing between the HRE velocity KDF ///(r, u, i) and the IKT 
approach. For this purpose we first remak that the representation of the HRE velocity KDF is actually non-unique. 
The proof of the statement is immediate. In fact, let us notice - preliminarily - that INSE can be equivalently 
represented in terms of the stochastic fluid flelds: 

{Zi(r,i,AV)}-{Vi(r,t)=V(r,f)+AV,pi(r,i),5T}, (120) 

where 

{Zi(r,i)} = {V(r,i),Pi(r,i),5T} (121) 

and AV = (AVi, AV2, AV3) E U = M.^ denote an arbitrary particular solution of the related initial-boundary value 
problems [i.e., Eas. (|137l) - (|139l) with the constant-entropy equation (J146I) ] and an arbitrary stochastic vector indepen- 
dent of (r, i) respectively. In particular, the NS equation for Vi(r, i) is manifestly 

^Vi(r,t)-fVi(r,i)-VVi(r,i) = 

= ^V(r, t) + (V(r, t) + AV) • VV(r, t) = F^ + AF^ , (122) 

with AFh being deflned as the stochastic vector fleld AFh = AV- VV(r, t). As a result, for an incompressible fluid the 
HRE velocity KDF, corresponding to Eq. ()122|) becomes ///(r, u— AV, t) = po5 (u — AV), which proves the statement. 
In view of these considerations let us now introduce the stochastic model defined by the set {AV S f/,(7(r,i, AV)}, 
with g{v,t, AV) being a suitable stochastic PDF [see Appendix B]. Due to its arbitrariness, it can always be identified 
with 

5(r,AV,0 = /(r,AV,i), (123) 

with / being defined by Eq.(fT3|) in terms of the IKT-statistical model {/, F}. This means that the corresponding 
KDF / is a particular solution of Eg. p^ obtained by prescribing the mean- field F in accordance with THM.l. As a 
result the following identity holds: 

(/H(r,u-AV,t))^v = /(i-,u,i), (124) 

where (•)^v i^ ^^^ stochastic average defined in Appendix B [see Eq. (|172|) ] with the stochastic vector a being identified 
with AV. In particular this implies, thanks to the correspondence principle (|31]) . that the variance of the stochastic 
velocity AV is prescribed as 



1 aaa2\ 

AV 



-AVn =Fi(r,i)- (125) 

"J / 1 
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As a further implication, when Ea. (jl24p is evaluated in the subspace of TTPs Ui, in view of Ea. (|93| and THM.3, it 
requires 

(/ff (r, uth - AV, t))^v = /(r, ^th, t) = -fi{r, Uth,t). (126) 

For all (r, Ut/i,i) G Fi x /, this yields also the relationship between the HRE velocity KDF and the conditional KDF 
fi(r,Uth,t) characterizing the statistical model of TTPs {/i,Fi}. 

A preliminary summary is in order. In the case of a NS fluid obeying the INSE problem, the following conclusions 
are reached: 

• In view of the constraint (I125p . g{r,AY,t) can be interpreted as the stochastic PDF taking into account the 
stochastic "thermal" motion of ITPs produced in a compressible thermal fluid by the kinetic pressure pi(r,i). 

• Eas. (|124l) and (|126l) permit a comparison between the HRE, IKT and TTP statistical models, {/i/,F}, {/, F} 
and {/i,Fi} respectively. We remark that such a comparison is always possible (and hence it holds also in the 
case of INSE). In particular, Ea. (|124p determines the relationship between the KDFs fH and / prescribed by the 
statistical model {/, F}. In addition, Eq. (|126p yields the analogous relationship with the KDF /i characterizing 
the TTPs statistics. 

• Remarkably, both conclusions follow by invoking a single suitable stochastic model. In both cases, in fact, the 
stochastic-averaging operator [see Eq. (|172l) ] which enters the l.h.s. of Eas. (|124p and (|126p is defined with respect 
to the same stochastic probability density ^(r, AV,i) prescribed according to Eq. (|123l) . 



• 



Eas. (|123p and (|124p - or equivalent (|126p - do not imply any restriction on the flow dynamics, i.e., on the (strong) 
solutions of the INSE problem. 

• Denoting {{))q the configuration-space average {{))^ = —r^Jdr, Ea. (|126p uniquely prescribes also the rela- 

o 
tionship between the corresponding spatial averages, i.e., 

(((/^(r, u,;, - AV, i))Av))o = ((/(r, u*., i)))o = I ((/i(r, u*,, t)))^ . (127) 

Hence, assuming that the initial frequency (((///(r, Uj/i — AV,io))Av))si ^^ prescribed, for the initial velocity 
PDF fi{r,Uth,to) the average {{fi{r,Uth,to)))^ remains uniquely determined too. 

Let us now address in detail the issues of the comparison between the HRE, IKT and TTP statistical approaches 
according to Eqs. (|124p and (I126P and the related closure condition problem arising in the HRE approach. 

For this purpose it is worth recalling that the aim of the HRE approach is actually to predict the time evolution, 
in the presence of turbulence, of suitable ensemble-averages of the KDF fH and of the NS fluid fields {V,p}, i.e., 
respectively (.fg(r u t)) and (V(r,t)), (p(r,t)). Here the brackets (•) denote a suitable ensemble-averaging operator 
(see for example [5^). In the case of so-called homogeneous, isotropic and stationary turbulence (HIST), this is 

required to commute with the differential and integral operators {Q} = < ^, ^, g^.a,. , ^j I '^^^I '^'^ (• -^^ pointed 

out above, in the context of the statistical description of turbulence, the operator (■) may be equivalently intended 



as a mean value in the probabilistic sense 57[, namely (•) = (•)^. Here (•)^ denotes again the stochastic average 
of the form (|172p [see Appendix B], prescribed in terms of a suitable stochastic probability density defined on the 
space Va of the stochastic parameters a. Thus, in this context the fluid fields are considered as stochastic functions 
dependent on a. In the case of INSE this requires letting {Zi{r,t,a)} = {V(r, i, a), pi(r,i, a), S'T(a)}. In particular, 
in validity of HIST this implies that necessarily the stochastic PDF g{r, t, a) must be taken of the form g = g{a). The 
corresponding statistical evolution equation for (/ij) is well-known and has been investigated by several authors (see 
for example [48j). In the case of the unforced NS equation, its explicit solution involves the construction of an infinite 
set of continuous many-point PDFs, coupled via the fiuid pressure, which obey a hierarchy of statistical equations, 
the so-called ML (Monin-Lundgren |5l|,|54]) hierarchy. 

The search of possible "closure conditions" for the ML hierarchy {Closure Problem for the statistical description of 
HT) remains - to date - one of the outstanding unsolved theoretical problems in fluid dynamics. Its solution involves 
in principle the search of possible alternative statistical models with the following features: 

Requirement ^1: it should be characterized by a finite number of (multi-point) PDFs. 

Requirement ^2: it should be determined in such a way that the complete set of fluid fields can be uniquely 
represented in terms of the same PDFs. 
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Requirement #3: the time evolution of the said PDFs is solely determined by a finite number velocity moments of 
the same PDFs {closure conditions). 

It is immediate to show that a possible candidate satisfying all of these features is provided by the TTP statistical 
model {/i, Pi}. In fact, on the basis of the theory developed above (see in particular THMs. 1 and 3), we conclude 
that, in the case of INSE: 

• The statistical set {/i,Pi} is realized in terms of the stochastic 1-point conditional PDF fi{r,Uth,t,a). 



• 



The fluid fields {Zi{r,t,a)} are all uniquely prescribed in terms of the same stochastic PDF /i(r, Uf/i,i, a) 
(correspondence principle). 

The time-evolution of the stochastic PDF /i(r, Uth,t, a) is uniquely prescribed by means of Liouville statistical 
equation [see the inverse kinetic equation (|94p]. 



• 



By assumption, such a kinetic equation depends functionally [see Axiom ^3 - Kinetic closure conditions] only 
on a finite number of velocity moments of the same PDF. 



Let us analyze how, in practice, the previous conclusions can be implemented in order to avoid the closure problem. 

In the customary HRE approach one is faced with the formidable issue of prescribing the multi-point PDFs which 
enter the Monin-Lundgren hierarchy. It is immediate to recognize that this problem arises specifically because of the 
treatment adopted for the fluid pressure in such approaches. In fact, because the pressure is not represented by a PDF 
velocity moment, but rather its contribution enters by means of a Green-function convolution integral, it follows that 
non-local (i.e., multi-point) contributions are necessarily introduced. The precise prescription of such contributions, 
however, remains undetermined, giving rise to the closure problem. Although several attempts have been suggested 
(see for example Ref.[53|), no definite solution exists to date. 

In contrast, based either on the IKT or TTP statistical descriptions, the problem can be given a consistent solution. 
This is reached as follows: 

1) By replacing the ensemble-averaged HRE-KDF (///(r, v,t)) = (//f (r, v,i,a))^ either with the corresponding 
IKT or TTP KDFs, namely (/(r,u,t)) = (/(r,u, t, a)}„ or (/i(r,Ut,„ t)) = {fi{r,Uth,t,a))^. 

2) In both cases the kinetic pressure is determined as a velocity moment of the corresponding KDF. This permits 
to overcome the closure problem arising in the HRE approach. In fact, non-local contributions due to multi-point 
PDFs do not appear anymore. 

3) By requiring that the KDF f{r,u,t,a), and respectively fi{r,Utfi,t,a), satisfy the Liouville equations (|14p and 
(1941) . These equations are determined imposing the requirement that they depend functionally only on a finite set 
of velocity moments of the same KDFs. Their numerical solution involves, at most, the determination of the local 
gradients of the corresponding fiuid fields. For comparison, the HRE-KDF obeys, instead, a Fokker-Planck statistical 
equation which contains non-local contributions due to the fiuid pressure which depend by higher-order multi-point 
PDFs (see for example Ref.[50]). 

4) Unlike the HRE approach, both IKT and TTP statistical approaches allow the stochastic fiuid pressure to be 
uniquely determined as a velocity moment of the relevant KDFs. Its stochastic average and fluctuating part follow 
simply by applying the stochastic average operator on the resulting expression. 

5) The choice of the stochastic variables a remains in principle arbitrary. Thus, they can be identifled in accordance 
with the speciflc stochastic model adopted (e.g., stochastic initial conditions, stochastic boundary conditions or 
stochastic volume force). 

6) Let us consider the comparison between the customary HRE statistical evolution equation (see again for example 
Ref.[50|) and the corresponding stochastic-averaged Liouville equations following from Egs. p^ and ([M]) upon applying 
the averaging operator (•)^. As pointed out in Section 6 and in Ref.'54], the latter equation can be approximated, 
locally in velocity space, in terms of a Fokker-Planck equation. The corresponding Kramers-Moyal coefficients however 
are different from those entering in the HRE equation. The remarkable features of our equation is, first, that unlike the 
HRE one, it recovers exactly the correct velocity moment equations obtained for the weight functions G = (1, v) (see 
Ref.[5J| on this issue) which also follow from the corresponding stochastic-averaged Liouville equation (Eqs. ([Ti| ) and 
(IMl)). Second, the Kramers-Moyal coefficients depend explicitly also on the stochastic pressure fiuctuations Sp, while 
the strict positivity of the KDF is warranted by the weak H-theorem (|113p following from the exact stochastic-averaged 
Liouville equations ([T4| or (p4|) . 

7) Another remarkable difference with respect to the HRE statistical equation is that the Liouville equations as 
well as the corresponding Fokker-Planck approximations presented here hold both for Gaussian and suitably-smooth 
non-Gaussian KDFs. Therefore, both IKT and TTP statistical models appear suitable to describe the non-Gaussian 
behavior arising in the statistical description of turbulence. 
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VIII. CONCLUSIONS 

A fundamental issue for Navier-Stokes fluids, is their characterization in terms of the dynamics of ideal tracer 
particles (ITPs) , and in particular of the sub-set of thermal tracer particles (TTPs) . Based on the formulation of an 
inverse kinetic theory for compressible/incompressible NS thermofluids, in this paper properties of TTP dynamics and 
a mathematical model for their description have been investigated. It is found that TTP dynamics can be uniquely 
determined both for incompressible (isothermal or non-isothermal) and compressible NS fluids described respectively 
by the INSE, INSFE and CNSFE problems. In addition it has been proved that the discovery of TTP dynamics 
allows for the construction of a reduced-dimension statistical model, to be identified with the TTP-statistical model. 
The latter is defined by the set {/i,ri}, in terms of which the self-consistent time evolution of the fluid fields is 
determined. 

Basic consequences of the theory mentioned in this paper concern the treatment of stochastic fluid fields arising in 
compressible/incompressible NS fluids by means of the TTP-statistical model. Here we have pointed out in particular: 

1. The formulation of TTP dynamics for the stochastic CNSFE problem, represented in two possible forma. The 
first one is provided by Langevin equations, which describe the dynamics of TTPs in fluctuating compress- 
ible/compressible NS fluids. These provide a mathematical model of tracer-particle motion in stochastic fluids. 
The second one is given by the corresponding Fokker-Planck description. 



2. The comparison with the statistical treatment of turbulence due to Hopf, Rosen and Edwards [2J-|26|, which 
holds in the case of incompressible isothermal NS fluids, has been carried out. As a result, the relationship 
between the HRE, IKT and TTP velocity PDFs has been pointed out. 

3. Finally, based on the statistical model {/i,ri}, a solution of the closure problem for the statistical description 
of HT has been proposed. 

Applications of the present theory are in principle several. They concern, in general, the dynamics of small 
particles (such as solid particles or droplets, commonly found in natural phenomena and industrial applications) 
in compressible/incompressible thermofiuids. The accurate description of particle dynamics, as they are pushed along 
erratic trajectories by fiuctuations of the fluid flelds, is essential, for example, in combustion processes, in the industrial 
production of nanoparticles as well as in atmospheric pollutant transport, cloud formation and air-quality monitoring 
of the atmosphere. 
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IX. APPENDIX A: DETERMINISTIC/STOCHASTIC NS THERMOFLUIDS 

Let us consider for deflniteness a viscous and generally non-isentropic thermofluid either incompressible or com- 
pressible, described in both cases by the fluid fields ^. In the following, we shall assume that the fluid flelds are 
deflned and suitably smooth in the existence domain fl x I, with fl and / denoting respectively an open subset of the 
3-dimensional Euclidean space R'^ and a subset of the real axis M. 



27 



A. A.l - Case of a compressible fluid: CNSFE 



Let us consider the case of a compressible thermofluid. Denoting by -^ ~ ^ + ^ ' ^ t^'-' fluid convective derivative, 
its fluid equations are identified with the so-called compressible Navier- Stokes- Fourier equations (CNSFE) 

-pV-V = 0, (f28) 

= Fff, (129) 

K, (130) 

(131) 



Dp 
Dt 

fv 

DT _ 
d 



p{r,t)^p{r,t,p,V,T), 



(132) 



where Eqs. (|128p - (|130[) denote the mass continuity, forced Navier-Stokes and Fourier equations respectively; further- 
more, the inequality (|13ip identifies the entropy law, customarily known as the 2nd principle of thermodynamics and 
Ea. (|132p is the equation of state for the fluid pressure. Here the notation is standard. Thus, in particular in Ea. (ll29p 
Fh denotes the NS acceleration 



P 



where the viscous stress tensor o_'{r,t) is assumed of the form 



^'(r,i) 



M 



VV+(VV)' 



-iv.v 



AlV • V, 



with /i, A > to be denoted as first and second viscosity coefficients. As a consequence 

2 



V-^'(r,i) = V- Af(vV+(VV) 



A 



;^^ 



v-v 



(133) 



(134) 



(135) 



In addition, f denotes the volume force density, which is assumed to be a suitable smooth vector field of the general 
form f — — V(/) -I- ffl, with (f>{r,t) and ffl(r, t; {Z(r, t)}) denoting respectively a scalar function (potential) and an 
additional non-conservative vector generally dependent of the fluid fields and in particular the temperature T(r,i). 
As an example, f can be identified with f = — pg [1 — kpTj+fm , where g and kp are real constants denoting respectively 
the local acceleration of gravity and the density thermal-dilatation coefficient, —pgkpT is the temperature-dependent 
gravitational force density and {jii{r,t;{Z{r,t)}) is a possible additional volume force density. Finally, the source 
term K in Eq. (|130p denotes the heat production rate, defined as 



Kir,t) 



1 



[Cp 



tP 



drk 



uJt — iPrn 



' Dt 



dr 



^(^ + ^-i%^) +A(V-V 




(136) 



where Cp is the heat capacity at constant pressure and a, Pt are suitable (dimensional) phenomenological parameters. 
In the case of an incompressible fluid, denoting po the constant mass density, v = p/ Po the kinematic viscosity and 
requiring X — a — /3t = 0, the previous equations reduce to the incompressible Navier- Stokes- Fourier equations 
{INSFE): 



p = Po> 

V • V = 0, 

D 

Dt 
DT 

iDt 

d_ 

di 



K, 
St >0, 



(137) 
(138) 

(139) 
(140) 
(141) 
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where now 



Fh = [Vp-f] 

Po 



i/V^V 



(142) 



K 



nJ^ + - f ^ + ^ - -5.,— 1 +A(V.V) +V-[kVT] 



(143) 



Furthermore the equation of state (|132p is replaced by the Poisson equation for the fluid pressure. This is obtained 
by taking the divergence of the NS equation (|139p . yielding 



V^p ^ -p\J ■ (V • VV) + V ■ f , 



(144) 



with p to be assumed non negative and bounded m Vt x I. Finally, we remark that Eas. (|137l) - (|14ip include, as a 
particular case, the treatment of isothermal fluids. This is obtained assuming an initial spatially uniform temperature 



T{v,to)=To, 



(145) 



requiring that for all (r,i) e fi x / the heat production rate K{Y^t) is identically zero in 57 x / and imposing, at the 
same time, the isentropic law 



P--' 



(146) 



Eas. (|137l) - (|14ip with the constraints (jl45p and (I146P are denoted as (isothermal and) incompressible NS equations 
(INSE). 



1. Equivalent forms of ttie Fourier equation 



To prove Eq. p30p with p36p let us start from the law of energy conservation equation. For a compressible viscous 
fluid this is HI: 






V ( ne + -pV^ +p] -Y -a' - kVT 



(147) 



with n{r,t) and p{r,t) denoting respectively the number and mass densities, while e{r,t),k > and T{r,t) > are 
the internal energy density, the thermal conductivity and the temperature. In terms of the convective derivative this 
delivers 



D 



n——e ■■ 
Dt 



-P- f-v^ 



V • [- Vp + V • (t' + fcVT] 
On the other hand from the NS equation [see Ea. (|129p ] it follows 



(148) 



P- (-v^ 

^Dt \2 



-V- Vp + V-f + V- V-CT'(r,i), 



so that Ea. (|148p recovers immediately the internal energy-transfer equation 

De 



Dt 



Se, 



1 



Se = JT + - l-pV • V + V • [kVT] + (VV) : CT'(r,i)| 



where 



(149) 



(150) 



Jt 



-Vf 



n 



(151) 
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and 

carry respectively the contributions to the due to the e-production rate generated by the vohime force f (i.e., external 
sources) and viscous energy dissipation. 

This equation, thanks to Eqs. (ll53p - (ll49p . can be equivalently cast into an equation for the temperature T of the 
form (I130p . In fact, for generally non-isothermal and compressible fluids e(r,t) can be taken such that 

D ( a \ DT / aT\ Dp 

Ea. (ll53p requires 

/ a \ DT / aT\ Dp 

" ['^ ~nP)-DT = -" [^^ - —) m + 

-p§-JlvA+V-[-Vp + V-a' + kVT]. (154) 

Introducing the definition (I15ip and invoking Ea. (|152p this yields 

a \ DT rr^^Dp 



{c,-%)^=nJT-{PTn-aT)^^-pV-Y- 



and hence Eqs. (|130p and (jl36p . If we introduce, instead, the representation of e(r,t) in terms of the local entropy 
s(r,i): 

D ^ Ds p Dlnn , ^ 

— e = T — + - , (156) 

Dt Dt n Dt ^ ^ 

Eq. ()148p delivers the local entropy equation: 

nT— = Ss, (157) 

Dt ^ ' 

S,(M).nJ. + V.[Wr| + f(f + f _1.,|1)^A,V.V)^ (158, 

with Ss denoting the local entropy production rate. 

2. Entropy law - Externally heated thermofluid 
Denoting the global thermodynamic entropy Srit) as 

Srit) = f drns, (159) 

n 

it follows identically 

;STit)^Jdr^^ins)^ Jdr 
n n 

I dv [-nV • Vs + 5 - sV • (nV)] = j dvS,. (160) 



di' 



' d d ' 

at at 



Hence, subject to the requirement: 
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dr^Jr > 0, (161) 



n 
usually referred to as (condition of) externally heated thermofluid - the entropy law (I13ip follows, with 



llSTit)^Jdr 



I kiyrj^ ii (dv dVk 2 dvA\\ . 

-Jr + —^;^ + -l— + — ~-S.,— +-(V.V) 



n 
denoting the global thermodynamic entropy production rate. 



(162) 



B. A. 2 - Deterministic and stochastic initial-boundary value problems 

The fluid equations defined by CNSFE are required to satisfy initial-boundary value problems defined by appropriate 
initial and Dirichlet boundary conditions 

Z{r,to) = Zo{y), (163) 

In particular, denoting / C R a suitable subset of the real axis, we shall require that Sxit) is defined and smooth for 
all i G / and that a smooth (strong) solution exists for the previous initial-boundary value problem in F x / (existence 
domain) . 

Finally, we shall assume that the fiuid fields {Z}, together with the volume force density f and the initial and 
boundary fields Zq and Z„, are all stochastic functions of the form (see Appendix B) 

Z = Z(r,i,a), 

Zw — Zyj [r,t,a)\gQ , 

with a G Vq stochastic variables independent of (r, t). Depending whether the previous functions (I164p are considered 
deterministic or stochastic, the previous initial-boundary-value problems defined either by: 

1. Eas. U28]) - il31\) with the initial-boundary conditions U63\) . 
or: 

2. Eqs.p^T^-pJ^ with HJ^ and the initial-boundary conditions il63\) . 

3. Eqs.[W^-(T^^ with pj^ and together with the constraints HJ^, PT^ (^""^d the initial-boundary conditions 

will be denoted, respectively, as deterministic or stochastic CNSFE, INSFE and INS E problems. 

C. A. 3 - Equivalent stochastic fluid equations 

For a prescribed stochastic model {g, Va}, with g — g{v, t, a) being a stochastic probability density on Vq,, in terms 
of stochastic decomposition ([M)) it is immediate to obtain the equations for the average and stochastic fiuid fields 
{Z)^, SZ and the corresponding initial-boundary value problem. For example, in the case of INSE, requiring that 
(•)q, commutes with the nabla operator V, Laplacian V^ and partial time derivative ^ operators, the fiuid equations 
for (Z)^ and 6Z, to be referred to as stochastic incompressible NS equations, become respectively 

^v\ =(F^)„-l[V(p)„-(f)] + ^.v2(V)„, (165) 

V ■ (V)„ - , (166) 

^V-/^v\ = 5Yh - — [V5p- 5i] + vVHY, (167) 

Dt \Dt / ^ po 

V • (5V = 0. (168) 
In particular Eqs. (ll65p - (ll66p identify the so-called stochastic- averaged INSE. 
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X. APPENDIX B - STOCHASTIC/DETERMINISTIC VARIABLES AND STOCHASTIC MODELS 

Let (Sj'SjP) be a probability space; a measurable function a : S — > Va, where Va C M.^ and fc > 1, is called 
stochastic (or random) variable. 

A stochastic variable a is called continuous if it is endowed with a continuous stochastic model {a G Vq,, (7(r, ijO;)}, 
namely a real continuous function g, called PDF on the set Va, such that: 

1) g is measurable, non-negative and of the form 

g^g{r,t,-); (169) 

2) if A C V^ is an arbitrary Borelian subset of Va (written A G B{Va)), the integral 

Pa{A)= f dagir,t,a) (170) 

A 

exists and is the probability that a G A; in particular, since a &Va, g admits the normalization 

dag{r,t,a) = Fa(V„) = 1 . (171) 

The set function Pa : B{Va) —^ [0, 1] defined by (J170I) is a probability measure on Va- Consequently, if a function 
/ : Va — > Vf C M™ is measurable, / is a stochastic variable too. 

Then we define the stochastic-averaging operator {■)^{aee also Refs.[2, ISJ]) as 

{f)a = (/(y, •))„ ^ / dag{r, t, a)/(y, a), (172) 

for any Pa-integrable function /(y, •) : V^ — > R, where the vector y is some parameter. 

The ensemble {g{r,t,a),Va} is denoted as stochastic model. Examples of stochastic model are represented by 
statistical models. In such a case the stochastic variables a are identified with hidden variables, i.e., variables 
from which the fluid fields {Z} depend only implicitly. For example, in the IKT statistical models {/, F} and 
{/i,Fi} (see THMs.l and 3, in Sections 4 and 6) the stochastic variables and PDF are prescribed letting a = u and 
g{r,t,a) = /(r, u,t) or g{r,t,a) = /i(r, u,t), with /(r, u, i) and /i(r, u, i) being defined, respectively, by Eqs. (fT3)) 
and (|89l) . Furthermore, in the two cases the set Va is identified with the velocity spaces f/ or Ui. 

A classification of stochastic models can be given in terms of the defining PDF g as follows. 

Definition - Homogeneous, stationary, deterministic and stochastic PDF 

The PDF g is denoted: 

a) homogeneous if g is independent of r, namely g = g{t, a); 

b) isotropic if 5 is a function of the form g — g{\r\ ,t,a); 

c) stationary if g is independent of t, i.e., g — g{r, a); 

d) deterministic if g is a distribution on Va of the ioTmg{r,t,a) = 6^'^' {a — ao), with 6^'^' {a — Uo) denoting the 
fe-dimensional Dirac delta on the space Va] 

e) stochastic if g is an ordinary function on the space Va ■ 
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